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ABSTRACT 



We report on the variability of 443 flat spectrum, compact radio sources 
monitored using the VLA for 3 days in 4 epochs at ~ 4 month intervals at 5 
GHz as part of the Micro-Arcsecond Scintillation-Induced Variability (MASIV) 
survey. Over half of these sources exhibited 2-10% rms variations on timescales 
over 2 days. We analyzed the variations by two independent methods, and find 
that the rms variability amplitudes of the sources correlate with the emission 
measure in the ionized Interstellar Medium along their respective lines of sight. 
We thus link the variations with interstellar scintillation of components of these 
sources, with some (unknown) fraction of the total flux density contained within 
a compact region of angular diameter in the range 10-50/xas. We also find that 
the variations decrease for high mean flux density sources and, most importantly, 
for high redshift sources. The decrease in variability is probably due either to an 
increase in the apparent diameter of the source, or a decrease in the flux density 
of the compact fraction beyond z ~ 2. Here we present a statistical analysis of 
these results, and a future paper will the discuss the cosmological implications 
in detail. 

Subject headings: galaxies: active — ISM: structure — radio continuum — radi- 
ation mechanisms: nonthermal 
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Introduction 



The discovery of centimeter- wavelength Intra-Day Variability (IDV) or "flickering" in 
Active Galactic Nuclei (AGN) by iHeeschenl (119841 ) initially raised concerns that some AGN 
possess brightness temperatures over six orders of magnitude above the 10 12 K inverse Comp- 



ton limit for incoherent synchrotron emission (e.g. iQuirrenbach et al.lll989l ). However, con- 
siderable evidence has now accumulated to demonstrate that interstellar scintillation (ISS) 
in the turbulent, ionized interstellar medium (ISM) of our Gala xy is the principal mechanis m 



responsible for the IDV observed in AGN, as was proposed by lHeeschen fc Rickettl (119871 ). 



Two more recent observational techniques provide compelling evidence for the preva- 
lence of ISS. Time delays of 1-8 min are observed in the arrival times of the flux den- 
sity variations between telescopes on different co ntinents for the three intra-hour variable 



sources B0405-385, B1257-326 an d J1819+3845 (IJauncev et all l2000l ; iBignall et all 12006 



Dennett-Thorpe fc de Bruynll2002l ). The delay arises due to the finite time required for the 
stochastic fluctuations associated with the ISM to drift across the Earth. A second observa- 
tional signature of ISS relates to the modulation of IDV variability timescales with a period 
of exactly one year. This arises because the Earth's orbital motion about the Sun contributes 
to the effective velocity with which the interstellar scattering material moves relative to an 
Earth-bound observer; the variations are slow as the Earth moves parallel to the material 
and fast as it moves anti-parallel to it. Annual cycles in IDV variability timescales are re 



ported in at least seven sources (e.g. iDennett- 



Jauncey fc Macquartll200ll ; IBignall et al 



2003 



horpe & de Bruvnl 20031 : iRickett et al.ll2001 



Jauncey et al.ll2003l ). including several whose 



long variability timescales preclude detection of time delays in the scintillation pattern over 
intercontinental distances. For many lines of sight through the ISM the slowest variations 
are expected in September if the motion of the turbulent material is comparable to the local 
standard of rest (LSR). 

The recognition of ISS as the dominant cause of IDV has not entirely alleviated the 
brightness temperatures problems posed by these sources. A source must be small to scin- 
tillate; in the weak scintillation case most frequently observed at frequencies near 5 GHz 
( jWalkerlll998l ). the source angular size must be comparable to or smaller than the angular 
size of the first Fresnel zone, Bp = ^c/2vt[L. Here L is the distance to the scattering region, 
which we will refer to as the screen even though in some cases it may be better described 
as a slab extending from the Earth out to distance ~ L. 6p is typically tens of microarcsec- 
onds for screen distances of tens to hundreds of parsecs, which implies source components 
with angular sizes two to three orders of magnitude finer than the scales probed by VLBI. 
The long time-scale over which IDV has been observed in some sources suggests that such 
scintillating components can be relatively long-lived despite their small physical sizes. 
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This paper reports on the results of a Micro- Arcsecond Scintillation-Induced Variability 
(MASIV) survey for IDV in AGN. This year-long survey conducted observations of between 
500-700 AGN over each of four epochs of three or four days duration in 2002 and 2003 at 
4.9 GHz with the VLA. The aim of the survey was to provide a catalogue of at least 100 
AGN which vary on timescales of hours to days to provide the basis of detailed studies of the 
IDV AGN population drawn from a well-defined sample. A description of the observations 
in epochs 2, 3 and 4 is presented in $2] as a supplement to des criptions of the first epoch 
observations and MASIV source selection in lLovell et al.l (120031 ) (hereafter Paper 1). In §3] 
we describe how the time series of flux density for each source in each epoch was classified 
as variable or non-variable. In §4] we describe how we have quantified the amplitude and 
timescale for the variations from an analysis of the structure function combined from all 4 
epochs and apply corrections for noise and other sources of flux density error. The basic 
hypothesis of the paper that the variations are predominantly due to interstellar scintillation 
is presented and examined in $5j including the influence of parameters of the interstellar 
medium ( §5.11 emission measure and galactic latitude) and of the sources themselves ( §5.21 
mean flux density, spectral index). We now have redshifts for more than half of the sources 
and we present the dependence of the variability on redshift in §5.51 - a result which shows 
that the MASIV survey provides a new cosmological probe. In §5.61 we consider whether 
the variability is intermittent over the four epochs. Our data are listed in Table H] and our 
conclusions are presented in §|6j 



2. Observations and Data Reduction 

The VLA observations took place over four periods during January, May and September 
2002 and January 2003 (see Tabled]). All epochs were 72 hours in duration except Septem- 
ber 2002, which included an additional 24 hours. The additional time in this epoch was 
added as an attempt to detect the slower variation expected in September due to interstellar 
scintillation caused by material moving at a velocity comparable to the LSR. 

All observations took place during array reconfiguration and in each case the array 
was being moved to a more compact configuration (except for epoch #1). In each epoch 
the VLA was divided into five independent subarrays. In the first epoch each subarray 
observed a subset of the 710 source sample. Subar rays 1 through 4 o bserved the "core 



578 sources of our compact, flat spectrum sample (lLovell et al.l 120031 ). namely the weak 
(105 < S8.6GHz < 130 mJy) and strong (58.6GHz > 600 mJy) sub-samples, while the fifth 
subarray observed a sample of intermediate flux density sources in two regions of the sky. 
For the next three epochs all sources previously observed in subarrays 1 through 4 were 
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reobserved, thus providing full coverage for the core sample. However, subarray 5 was re- 
dedicated to observing a smaller number of sources comprising all objects found to be variable 
in subarray 5 during the first epoch, as well as the most rapid variables found during the 
first epoch which required faster sampling in order to determine time scales. In this paper 
we compare results for the core sample from subarrays 1 through 4 only. 

Shadowing was a serious concern at low elevations in the more compact VLA config- 
urations. The impact of shadowing was minimized by assigning 5 or 6 antennas to each 
subarray in such a way that at any given time a source could be observed with at least 3 
unshadowed antennas. Any data from baselines containing shadowed antennas were flagged. 
Antennas that had either not moved during reconfiguration or were moved but had pointing 
solutions already applied were assigned to subarrays 1 to 4 where possible. 

The core sample was divided into four roughly equal parts by Declination and a subarray 
assigned to each. The Declination ranges for the four subarrays were 0° < 5 < 14°, 14° < 
5 < 34.08°, 34.08° < 5 < 49°, 5 > 49°. The Declination boundary between the second 
and third subarray was set to the latitude of the VLA to avoid long slews in azimuth when 
changing between sources transiting north and south of the zenith. 

Each subarray was scheduled so that every source was observed for one minute every ~2 
h while it was above an elevation of 15°. We used the standard VLA frequency configuration 
for continuum 4.9 GHz observations (dual polarisation and two 50 MHz bandwidth IFs 
per polarisation) and a 3.3 s integration time. For flux density calibration, each subarray 
observed B1328+307 (3C286) and J2355+4950 every ~ 2 hours. B1328+307 is the primary 
flux density calibrator for the VLA and J2355+4950 is a GPS source, not likely to vary over 
short timescales, and is monitored regularly at the VLA as part of a calibrator monitoring 
program. 

Following the observations we calibrated the data in AIPS using the standard technique 
for continuum data. The task FILLM was used to load the data where corrections were made 
for known antenna gain variations as a function of elevation and for atmospheric opacity. 
Upon inspection of the data it was clear that there were residual time- dependent amplitude 
calibration errors. We ascribe this to the fact that in each epoch, some VLA antennas had 
recently been moved and their pointing calibration observations were not complete. The 
residual pointing errors may depend on azimuth and elevation so we chose several bright, 
non- variable sources in each subarray at a range of RA as gain calibrators for surrounding 
sources. Precautions were taken to ensure that the calibrators themselves were not variable: 
if a given calibrator caused the majority of sources against which it was applied to vary, 
then another calibrator was chosen. These sources, typically recommended VLA calibrators, 
were drawn from our source sample. The numbers of sources used as secondary calibrators 
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in the four epochs were 42, 33, 20 and 36. On the epochs that a particular source was used 
as a calibrator, it is by definition non- variable and was excluded from the structure function 
analysis discussed below. 

Following calibration, the data for each source were inspected and occasional outlying 
samples were flagged. The data were then incoherently averaged on a one-minute timescale 
over all baselines. For all sources, the formal errors obtained were less than those estimated 
due to the residual constant and fractional errors discussed later in this paper. Incoherent 
averaging was chosen because, on the assumption that all sources are unresolved, the phase 
should be zero and any residual phase errors in the data would artificially reduce the average 
flux density. In low signal-to-noise regimes an incoherent average can induce an upward bias, 
as visibility amplitudes foll ow a Rice distribution, which has a mean that is systematically 



higher than the true mean ([Thompson et al.ll200ll ). In the case of our observations though 



the signal-to-noise ratio of our observations is sufficiently high that this effect is negligible 

As the VLA array configurations became more compact, our data became more sensitive 
to extended structure in the source and in nearby objects. The changing response of the 
VLA with time due to this structure can appear as variability in a light curve. Fortunately, 
as our observations were scheduled in sidereal time and were repeated at the same times 
every day, false variability due to structure appears as a repeating pattern with a period of 
one day. In general this was easy to recognise but in some cases imaging was carried out to 
verify contaminating structure. For the purposes of our analysis, such sources were removed 
from our sample. A total of 102 sources were removed due to structure or confusion, and 
one was removed due to an error in our initial sample selection process. Following these 
removals, we are left with a sample of 475 point sources common to all four epochs. 



3. Classification of Variables 

We used two separate approaches to analyze the flux density variations. This section 
describes how we classified the sources based on the apparent modulation index of their 
intensity. The following section describes a structure function analysis. The virtues of these 
two techniques are complementary: the first is based on a conservative but robust criterion 
of source variability, while the second uses a simple statistical estimator which allowed us to 
quantify the amplitude and time scale of the variations. 

In order to ascertain whether an individual source is variable it is necessary to under- 
stand the sources of error inherent to the measurements. There are uncertainties in the 
individual measurements due to calibration errors causing a fractional error, p, and additive 
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errors, s, due to thermal noise and confusion. Calibration errors include contributions from 
antenna pointing errors, system gain variation between the observations of flux calibrators 
and variable atmospheric absorption. Since they are a small percentage of the mean source 
flux density, S, they can be approximated as additive and added in quadrature to the noise 
as given by 



Here cr er r,s,p is the rms error in each flux density estimate normalized by the mean for each 
epoch. In our initial analysis (Paper I) we estimated s = 1.5 mJy and p = 0.01. 

An initial inspection classified sources based on a their modulation index, defined as 
the rms of the three days observations divided by the mean flux density, computed for 
each epoch. A source was identified as variable if its modulation index exceeded twice the 
expected contribution from the measurement errors, 2er err;SjP , as in a x 2 test. However, 
direct inspection of the data revealed that many of the slower variables, that is sources 
with variability time-scales longer than three days, were not detected as variable using this 
criterion. 

We therefore introduced an alternative variability criterion, which classified a source as 
variable if the modulation index of its daily average flux densities exceeded 2<7 err . This process 
yielded more detections, but visual inspection again revealed that many sources that clearly 
exhibited variability were undetected by either test. In particular, examination revealed 
that our two criteria do not detect variability in those sources with low-level monotonic flux 
density changes during the 3-4 d duration of our observations. This is because the \ 2 statistic 
used in our two selection criteria is not an ordered statistic and is therefore sub-optimal in 
detecting such a low level trend. 

Our selection criteria were therefore augmented with a visual inspection of the remaining 
light curves. Inspection was performed by two of us independently and the results compared. 
Each source was considered non-variable unless otherwise demonstrated, and any source 
where there was disagreement on its classification was reviewed; we adopted the conservative 
approach of classifying as non- variable any source is which no final agreement was reached on 
the classification. In the left panel of Figure [T] we show a scatter plot of the raw modulation 
index m raw of all sources against their mean flux density from epoch 1, using differing symbols 
for variables and non-variables. The predicted error in a single flux density measurement 
from equation ([1]) is shown by the dashed line, which roughly separates the variable from 
non-variable classifications. This emphasizes that our classification is determined largely by 
m raw . Table H] lists all the raw modulation indices and variability classifications of all the 
sources and Figures [2] and [3] show some sample light-curves from the survey. 
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In Table [2] we list the number of sources classified as variable in zero, one, two, three or 
all four observing epochs and associated percentages. In this table and in all of what follows 
we define a set of 443 sources from the 475 sources that excludes those used as secondary 
calibrators in two or more epochs. Note that 12% of the sources were seen to vary during all 
four epochs. With any analysis of a large number of observations, however, false positives, 
i.e. sources that are incorrectly classified as variable, are a significant concern. Since the 
visual classification was very close to the two sigma criterion, our variability classification 
is reliable with ~ 95% confidence. Thus we give in columns 3 and 4 the predicted fraction 
of sources misclassified based on the null hypothesis that all 443 sources are non-variable. 
Evidently these are negligibly small relative to the observed fractions of variables in two or 
more epochs. 

The large numbers of sources classified as variable on multiple epochs firmly establishes 
that our classification process is reliable. If all 443 sources were non- variable, then the number 
expected to be classified variable on 2, 3 or 4 epochs is a mere 6.7 or 1.4%, compared with 
the observed 192 (43%). Even if our classification was 90% reliable, then the number of false 
positives with 2 or more detections remains less than 5%. 

The fraction of sources that can be reliably classified as variable can be deduced using 
Table [2j Denoting N as the actual number of non- variable sources, T the actual number 
of one-time variables and F the number of non- variables misclassified as one-time variables, 
then the apparent number of non- variables 161 = N — F and from Table [21 F/N = 0.172. 
This yields iV = 194, F = 33 and hence T = 90 — F = 56, and the corrected total number 
of variables is 192+56=248. Therefore the fraction of sources that exhibited variability in 
one or more epochs in our surv ey is 56%. This value is significantly higher than the 15- 20% 



found in previous IDV surveys (IQuirrenbach et al.lll992l ; iKedziora-Chudczer et al.ll2001al ) . In 



comparing with other surveys one must be careful to consider the selection criteria applied. 
As described in §2] we started with 710 sources selected from spectral index and mean flux 
density criteria which was reduced to 443, when we excluded those used as calibrators in 
more than one epoch and those exhibiting resolution effects, raising the percentage of variable 
sources. In addition the large number of one-, two- and three-time variable sources raises 
the percentage. We believe that this is due to intermittency in the IDV phenomenon, which 
we describe by a simple model in §5.61 

For the purposes of subsequent analysis we conservatively define as "non- variable" those 
161 sources that showed no variability in any of the four epochs, and we define as "variable" 
those 192 sources that showed variability on two or more of the four epochs. With these 
definitions we have two large and reliable samples each of approximately 200 sources, where 
the non-variables act as a control sample for the variables. Each was drawn from the same 



- 9- 



selection criteria and cover the same overall area of sky. 



4. 



Structure Function of the Variations 



Here we discuss how we quantify the flux density variation for each source using the 
structure function (SF) of each time series, 



where Sj is a flux density measurement normalized by the mean flux density of the source 
over all four epochs and N T is the number of pairs of flux densities with a time lag r binned 
in 2-hour increments. The SF is a statistically reliable estimator which can be modeled even 
for short data spans. It is defined independent of any variability classification. Examples 
are shown in the lower panels of Figure [2] and [31 The SF is preferable to the autocorrelation 
function for short data spans, which can be badly biased by a poor estimation of the mean. 

For an idealized observation of stationary stochastic variations spanning a time much 
longer than their characteristic time r c h ar? D{r) rises with time lag and tends to saturation 
at twice the true variance. However for our observations r c har is typically more than two 
days and saturation is rarely seen. We thus chose D(r = 2d) as a standard characterization 
of the intra-day variations because r = 2 d is the maximum lag out to which our structure 
functions contain reliable information for the 72 hour observing sequences. Note that we 
estimate D(t = 2d) by combining data from all four epochs, which makes it a more stable 
statistic under the hypothesis of stationary stochastic variations. Hence we do not use this 
statistic to determine whether a source's variability differs over the four epochs. 



Additive flux density errors that are independent of the ISS (or intrinsic) variations 
contribute an additive term to the SF. If the errors are independent from one sample to 
the next (i.e. "white"), the error contribution to SF would have a mean value -D noise = 2<7g rr 
independent of time lag. While random pointing errors will be white, there may be non-white 
errors caused by systematic pointing errors due to recent antenna relocations and residuals 
of gain and atmospheric variation, which may contribute a term that is a function of time 
lag. Nevertheless we now proceed by assuming that all the errors are white and thus the 
errors add a constant (-Dnoise) (plus random variations about the constant due to estimation 
error) and postpone discussion of possible non-white measurement errors. 




(2) 



4.1. 



Measurement errors 
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Our goal is to subtract D noise from the raw SF D(t) in order to estimate the SF of any 
true variations in a source's flux density. But first we must consider how best to estimate 
-Dnoise- One estimate is 2a^ rr from equation (jTJ), which depends on the constants s and p. 
Another estimate of D noise comes from the SF itself - evaluated at its shortest time lag 
(2h). This would be an unbiased estimate for -D no i S e if all real flux density variation had 
time scales much longer than 2 h. Hence we examined single-epoch estimates of D(t = 2h), 
which we plot as an equivalent modulation index m 2 h = ^/0.5D(t = 2h) in the right panel 
of Figure [U The sources classified as variable include a substantial number with large values 
of m,2h, which are due to rapid flux density variations stronger than expected for noise. Note 
in particular the highest point which is quasar J1819 + 3845 which shows large amplitude 
rapid ISS. Its timescale is typically shorter than our 2 h sampling and so its SF has already 
saturated at 2 h and its variations are "white" in our sampling. Similarly the other sources 
with elevated m2h are probably due to ISS with short time scales. 

Now consider the non-variables, which are plotted as pluses and provide a set of sources 
with low or zero variation in epoch 1, which are useful in studying the noise processes. 
In the right panel of Figure Q] the mean of the plus symbols lies significantly below the 
dashed line (equation (JTJ with s = 0.0015 Jy and p = 0.01) particularly for the higher flux 
density sources. The solid line with s = 0.0013 Jy and p = 0.005 provides a better model 
for the noise as discussed in §4.21 In the absence of any real variations the estimates m 2n 
should be scattered equally above and below their mean value. Comparing the right and left 
hand panels we see that m raw is typically higher than m 2 h even for the non-variables, which 
suggests that m raw is increased due to low level variations with a timescale longer than 2 h. 

The discussion above reduces the choice for estimating -D no i se to either D(2 h) or the 
value corresponding to the solid line D exXiSp = 2a^ rr with p = 0.005. Whereas it is appealing 
to use the observable D(2h), any rapid real flux density variations can contribute to D(2 h) 
which thus overestimates Aioise- Thus we adopted D err s>p as an estimator for -D no i se , which 
states that the measurement errors are well described by equation([T|) but with calibration 
errors contributing at a 0.5% rather than a 1% level. However we discuss a slight revision 
of this in the next section. 



4.2. Structure Function Correction and Fitting 

For each source at each epoch we computed the raw SF as defined in equation (j2J). As 
the time lag increases the number of available lagged products (N T ) drops and so increases 
the error in the SF. A threshold for plotting an estimate for D{t) was set at N T > 20% of 
the total number of data samples in that epoch. As the example in Figure [2] shows, this 
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gives estimates clustered at lags 2-8 h and near 24 h and near 48 h. In order to characterize 
the variability amplitudes we initially estimate the SF at a lag of 2 d from the mean SF and 
its estimated error, calculated using the values at lags in the range r = 48 ± 2 hours. The 
overplotted model is described in §4.31 

First consider the equivalent 2-day modulation index ni2d = a/0.5-D(2 d) (without noise 
correction) plotted against the mean source flux density S for epoch 1 in Figure HI It should 
be compared with the left and right panels of Figure [T] for the non-variable sources. The 
comparison makes it clear that rri2d > ^2h for almost all of the no n- variable sources. Indeed 
there is a close correspondence between rri2d and m raw . Thus it is clear that there are low level 
flux density variations on a timescale longer than 2 h in the sources classified as non-variable. 
In order to investigate what causes these we examined the difference AD = D(2 d) — D(2 h) 
for the non-variable sources from each epoch. 

We examined how AD depends on mean flux density, Galactic latitude and Ha emis- 
sion - quantities which we find in §|5] influence _D(2d) for the variable classifications. The 
results showed that the mean of AD in each epoch was significantly higher for the sources 
weaker than 0.4 Jy than for those stronger than 0.4 Jy. However, it showed no significant 
dependence on Galactic latitude or Ha emission. Hence the process responsible for the low- 
level variability in the sources classified as non- variable is not ISS and is unlikely to have an 
astrophysical origin. 

We consider the most likely cause to be confusion, which can be due either to extended 
source structure which is partially resolved on the baselines of each sub-array or to low 
level confusing sources in the primary beam. We had already eliminated obvious cases of 
confusion by removing sources whose light curves showed clear daily patterns of variability 
which repeated in each day of a 3 day sequence. Since our SF estimation is from the light 
curves normalized by S an increase in AD at low S is consistent with the effect of low level 
confusion characterized by a certain rms in Jy. This is supported by the VLA documentation 
that gives 2.3 mJy as the brightest source expected in a single antenna beam at 5 GHz. 

We examined the shape of the mean D{r) for high and low 5* sources in each epoch. 
Evidence for the effects of confusion was found in the consistent minima in D(t) near time 
lags t = Id and 2d. However, there was also a component rising from the noise floor at 2h 
and typically saturating between 12 and 24h. All three components were substantially larger 
for the weak group of sources than for the strong sources. When averaged over four epochs 
the noise floor was, respectively, at ~ 3.0 x 10~ 4 and ~ 0.61 x 10~ 4 and the averages for 
AD were, respectively, ~ 2.8 x 10~ 4 and ~ 1.1 x 10~ 4 for these groups (mean flux density 
0.13 Jy and 1.4 Jy). By equating these average SF amplitudes to 2o~\ rT using equation ([I]) 
we estimated s = 0.0013 Jy and p = 0.0073. The s values for each epoch ranged from 
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0.001-0.002 Jy and for p 0.003-0.01, but no consistent patterns were seen versus epoch. 

We also considered the effect of long-term variations as characterized by the variation 
between epochs of the mean flux density from each epoch. These slower variations are 
mostly intrinsic and might contribute to a trend within the 3 or 4 days of each epoch. So 
we calculated the rate of change of flux density from the average of the magnitude of the 
differences in flux density between neighbouring epochs, divided by the number of days 
between epochs and the resulting magnitude of D(2d) due to such a trend. Surprisingly, 
for all sources this was smaller than the noise-corrected D(2d) except for those which were 
negative as a result of noise subtraction. Further the highest of them was 0.0001, which is 
one quarter of the threshold value. So we conclude that long-term variations did not make 
a significant contribution to the variations observed within the epochs. We note, however, 
that we can usefully estimate the structure function on time lags of 3 and 6 months from the 
MASIV survey data, which may be useful in studies of intrinsic variability. Thus we include 
the epoch averaged flux density for each epoch in the data table HI 

The foregoing studies reveal that the apparently non-variable light curves include not 
only white noise but also a low level contaminating process whose rms values can be ap- 
proximately characterized by equation ([1]). While we can quantify the white noise process 
by this equation, the non-white contamination is not well enough understood to be reliably 
characterized by an SF which could be subtracted from the estimates for all sources at each 
epoch. So in order to reduce the effect of these non-white variations, we increased the esti- 
mated .Dnoise to to be subtracted by using the values s = 0.0013 Jy and p = 0.007, and we 
also set a threshold on D(2 d) = 4 x 10~ 4 below which the noise-corrected SF values may be 
contaminated and therefore should only be interpreted as upper limits to the true SF of flux 
density variations. 

It is of interest that, since confusion effects will be precisely repeated at 24 h intervals, 
the samples of the structure function at 24 and 48 hr will be unaffected by confusion. The 
results we present in subsequent sections are from fitting the SF over all time lags, since the 
fit makes better use of the data. However, we also analyzed the single sample estimates of 
D(24h) and £)(48h) and found very similar results, though with somewhat worse statistical 
errors. As an extra precaution we re-reviewed all the light curves and SF plots for 24 h 
periodic patterns and found 34 sources that might be contaminated at a level near the 
threshold. For these sources we used the lower of D(2d) from the fit and that estimated 
from 48 ± 2 h. 
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4.3. Variability Timescales 

Although we estimated the SF as described above for each source at each epoch, the 
single-epoch D{2 d) is only based on about one independent sample of a 2 d vari ation and so it 



has a large statistical error - such that its rms error is about equal to its mean (IRickett et al. 



20001 ). Thus we do not attempt to evaluate the variability amplitude on a source by source 
basis for each of the epochs. Rather we average the SF for each source over all 4 epochs 
(after subtracting the -D no ise as defined in the previous section). Hence the SF results are 
insensitive to any intermittency or annual changes in ISS due to the effects of the Earth's 
velocity. 

We then fitted the following simple model to the SFs: 



T 



D(T)=2m 2 — , (3) 

T + T char 

where 2m 2 is the amplitude at which the function would saturate and T cria r is the character- 
istic timescale where the SF reaches half of its saturation. The motivation for this form is 
described in Appendix |A] It approximates the form expected from ISS caused by a turbulent 
interstellar medium uniformly distributed through a thick scattering region, as opposed to 
turbulence confined to a thin layer (see equation flA4j) ). We estimated two parameters only: 
the timescale r char and the value of D(t = 2 d) for each source. It should also be noted that 
a light curve that is dominated by a linear trend in flux density gives rise to a parabolic 
SF, which is not well fitted by equation (j3J). The value of D{r = 2d) in such a case will be 
somewhat underestimated. 

An example of fitted SF is shown in Figure [21 The points show D(t) increasing (noisily) 
with time lag r but not reaching saturation. Since the timescale is defined at half the 
saturation value it is poorly constrained in this example: r cria r — 1-0 ± 0.5 days. However, 
from the same fit D(t = 2d) = (44 ± 3) x 10 -4 , which is quite well constrained. With 
observations limited to 3 days (4 days for epoch 3) it is not possible to estimate timescales 
longer than about 3 days. However, in those cases it was possible to recognize that the 
characteristic timescale is longer than 3 days from the shape of the structure function. Two 
other examples are shown in Figure [3], in which there is evidence for faster variations. For 
source J0949+5819 the variations are very strong in epoch 1 and much weaker in epoch 3. 
From the epoch-average we find r char = 0.02 ± 0.05 days and D(t = 2d) = (19 ± 2) x 10~ 3 . 
For J1328+6221 the variations are more consistent over the epochs with r = 0.2 ± 0.2 days 
and D(t = 2 d) = (10 ± 1) x 10~ 3 . Given our 2 hr sampling and the typically large fractional 
errors in the timescale we have simply classified the timescale into fast r criar < 0.5 d, medium 
0.5 < r c h a r < 3d and slow 3d < T criar . We also looked for any correlation between the 
timescale and D(t — 2d) but found no consistent pattern. 
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During the visual examination of the light curves for each source at each epoch, the 
timescales were estimated by counting the number of inflection points (i.e. change in sign of 
the derivative) for those epochs classified as variable. Since in the visual examination there 
was an effective smoothing, inflection points due to noise-like deviations were not counted. 
The majority of sources were found to show none or at most one inflection point indicating 
variability timescales that are predominantly longer than 3 days. The observed distribution 
of inflection points is shown in Figure [51 Only a small number of sources (20%) showed 2 
or more inflection points. A comparison of the distribution of inflection points for the weak 
and strong sources revealed no significant difference between the two classes. Overall, the 
distribution of timescales was statistically the same for each epoch, remembering that epoch 
3 was four days rather than the three days of the other epochs. An important conclusion from 
the timescale study is our 3 or 4 d lightcurves commonly underestimate both the timescale 
and true modulation index for many of the sources. 

The annual cycle reported in a number of IDV sources is due to the changing relative 



veloc ities of the Earth and the ISM responsible for the scattering (IMacquart fc Jauncey 



20031 ). If the ISM velocities follow the Local Standard of Rest, many sources would be 
expected to exhibit slower variations in the third quarter of the year (i.e. during the third 
epoch), and hence may more easily be missed because of the lengthened time-scales. Figure [5] 
shows the numbers of variables found at each of the four epochs. A Chi-squared contingency 
test shows no evidence that the numbers differ from the mean in any epoch, even though 
epoch 3 lasted for four rather than three days. The uniformity of variable numbers in 
each epoch suggests a lack of evidence for a third-quarter slow-down, and it follows that the 
majority of the scattering material is not moving at the LSR. This is perhaps not unexpected; 
both PKS 1257-326 and J1819+3845, the two sources for which reliable screen velocities 
have been measured, have measured screen velocities that differ significantl y from the LSR 



(jBignall et al.ll2006l : iDennett-Thorpe & de Bruvnll2003l : iLinskv et al. 1120071 ). 



In summary we used the visual analysis to classify each source at each epoch as variable 
or not variable. We computed SFs for all sources and then examined the SFs of those 
classified as non-variable in order to quantify the measurement errors. We were able to 
correct the SFs by subtracting a constant versus time lag due to errors that are independent 
over the 2 h sampling and are characterized by equation [TJ In addition we found a low 
level contaminating process with a timescale of 1-2 d, which we suggest is due to low-level 
confusion with an rms of 1-2 mJy The SF of the slower contamination could not be reliably 
estimated and so sets a limit on the minimum detectable variation in flux density. By fitting 
a simple curve to the epoch- averaged and noise-corrected SFs, we estimated D(t = 2d) 
for each source. The contamination is minimized by requiring this quantity to be above a 
threshold of 4 x 10 -4 for a useable estimate of the timescale r c h ar - Most sources were classified 
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as slow variables. 



4.4. Comparison of SF with the visual variability classification 

We now compare the SF analysis with the variability classification from £J3] Figure [6] 
shows the number of epochs in which a source was classified as variable plotted against its 
value of .0(2 d) obtained as described from a fit to the structure function of the cumulative 
data from all 4 epochs. The large circle shows the mean values for each group of sources. 
While it is clear that the sources with higher .0(2 d) were classified as variable more fre- 
quently, there is a very wide distribution in the rms level of the variation over 2 d among 
the sources. The vertical line marks the threshold .0(2 d) = 4 x 10~ 4 above which we have 
made a timescale estimate. Values below this should be regarded as upper bounds in view 
of the possibility of low level confusion. 

Table [3] lists the source counts sorted by the number of "variable" epochs for the 443 
sources (as always excluding those used as calibrators in more than one epoch.) It also shows 
the mean values of D(2 d) and the numbers of sources above and below the SF threshold. In 
total 37% of them are above the threshold versus 45% from the variability classification on 2 
or more epochs. In the latter classification process we do not characterize the rms amplitude 
but attempt to quantify any intermittency in the phenomenon. In contrast the SF analysis 
quantifies the rms amplitude over 2 d averaged over all 4 epochs. 



5. Interpretation as Interstellar Scintillation 



We now examine our basic hypothesis, stated in the introduction, that the variations 
in flux density detected in the MASIV survey are caused by interstellar scintillation (ISS). 
The extremely small diameters of pulsars revealed the ISS phenomenon almost as soon as 
pulsars were discovered and they still provide the best information on the distribution of 
small scale structure in the ionized interstellar medium - indeed they provide a calibration 
of the ISS phenomenon. The pulsar observations have been c ombined into a m o del fo r the 
distributio n of electron dens i ty in the interstellar medium by iTaylor fc Corded (119931 ) and 
revised by ICordes and Laziol (120051 ) . 



Because the refractive index of an ionized medium varies with radio frequency, there 
is a transition frequency (/ w ) above which the scintillation of a point source, like a pulsar, 
is weak in the sense that its scintillation (modulation) index (m pt ) is less than one. This 
frequency is on the order of 5 GHz but depends on the strength of the turbulent fluctuations 
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in electron density on the given line of sight (IWalkerlll998l ; ICordes and Lazioll2005l ). Above 
/ w the ISS of a point source has a single timescale approximately given by tp = t f/V, where 
V is the effective velocity of the Earth through the ISS diffraction pattern and rp = \J L\/2 
and L is the typical distance to the scattering region. 



7T 



In our observations the angular diameters of the extra-galactic sou rces are considerably 
larger than those of pulsars and so their ISS is heavily quenched. See iRickett et al.l ( 120061 ) 
for a discussion of how the "low-wavenumber approximation" can be applied for quenched 
scintillation even below the transition frequency. If we approximate the scattering as taking 
place in a thin region of the Galactic plane, we obtain si mple exp r ession s for the reduction 



Rickettl (]1986D ). In Appendix 



in scintillation index and lengthening of timescale (e.g. 
we apply the same simple "screen" model to the structure function analysis and obtain 
expressions for how the observable _D(2d) might vary with angular size of each source and 
on distance to the scattering region and the level of turbulence on each line of sight. Of 
course the level of ISS also depends on properties of the source - in particular the fraction 
of its flux density in its most compact component and on the effective diameter of that 
component. 



5.1. Galactic Dependence of ISS 



We start by comparing our D(2d) results with the emission measure (column density 
of the square of the electron density) as estimated from observations of Ha emission. We 
find the intensit y of Ha emission (i n Rayleighs) from the WHAM Northern sky survey on 
a 1 degree grid (IHaffner et al.l I2003T) n e arest to each source. We use the intensity summed 
over all velocities, which IHaffner et al.l (120031 ) interpret as proportional to the ISM emission 
measure on that line of sight, assuming the temperature of the emitting gas does not vary 
by a large percentage. We expect the level of ISS to be rela ted to the emissio n measure on 
that l ine of sight, as described by lSpangler fc Corded (119981 ) and observed by lRickett et al. 
(l2006h . 



Figure [7| plots D(2d) against the WHAM Ha emission (in Rayleighs). Though the 
scatter plot in the top panel shows little obvious trend, the bin averages in the middle panel 
show a clear upward trend with emission measure, which establishes ISS as the dominant 
cause of the variability in the MASIV survey. We stress the complete independence of the 
two data sets in this figure and that the bin averages are independent of any threshold set 
on D(2d). We exclude the extreme IHV source J1819+3845 from this and subsequent bin 
average plots because its D(2d) value is 0.25 which is so much higher than the next highest 
at 0.015 that it distorts the mean and the variance within its bin. 
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The bottom panel shows that the fraction of slowly scintillating sources clearly increases 
with emission measure and vice-versa for the fast scintillators. This finding that the longer 
timescales occur when seen through greater column density of electrons is consistent with 
enhanced ISS from strongly ionized regions of the ISM, which are typically at low Galactic 
latitudes and at greater distances L. An increase in L increases the scale of the scintillation 
pattern which slows the scintillation time - see Appendix [A] 

Since the strength and effective distance of the scattering layer depends on Galactic 
latitude, we also expect a dependence of ISS on Galactic latitude. For our visual classification 
of sources (as variable or non-variable) we asked the simple question "are their latitude 
distributions the same?". A Chi-squared contingency test dividing the sources into two 
samples, a low latitude sample, |6| < 40 degrees, and high latitude sample, |6| > 40 degrees, 
shows that the two distributions differ at the 98% confidence level. There are fractionally 
more variables at low latitudes than there are at high latitudes, supporting ISS as the origin 
of the intra-day variability. 

For the structure function analysis we si mply plot D(2d) aga i nst th e Galactic latitude 



of each source, in a fashion similar to that of iHeeschen fc Rickettl (119871 ). The upper panel 



of figure M is a scatter plot, differentiated by the timescale group (fast, medium or slow). 
The middle panel averages D(2d) into 30 degree wide bins for all sources, which as already 
noted is independent of the threshold on D(2d). There is a low level of scintillation above 60 
degrees, increasing in the mid range (30-60 deg), in both northern and southern hemispheres. 
However in the low latitudes (0-30 deg) the ISS increases in the south of the plane but 
decreases north of th e plane. While the fi gure is in reasonable agreement with the latitude 
dependence found by Rickett et al. ( 20061 ) in their analysis of the modulation index of 146 



flat-spectrum sources observed at 2 GHz with the Green Bank Interferometer, there are 
competing effects in our MASIV survey since data were sampled for no more than 4 days. 

As the latitude b decreases both the distance L to and the path length through the 
scattering medium increase (oc cosec b). The increased path length makes the scintillation 
stronger at lower latitudes so that the modulation index should increase. However, the 
increased distance L increases the scale of the scintillation pattern which slows the scintilla- 
tion time so that the structure function will saturate at times longer than 4 days, causing a 
decrease in D(2d). The combination of these two effects requires careful modeling. 

The model for the structure function described in Appendix A is a starting point for 
analyzing the effects of Galactic latitude. Equation flA50 shows that as the distance L 
increases D(2d) should decrease. However, the equation assumes that the scintillation index 
of a point source m p t = 1, omitting any increase due to the longer scattered path length 
at lower latitudes. In a more realistic model m pt may be less than one looking out of the 
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Galactic plane (\b\ ~ 90 deg) and should increase with decreasing latitude, reaching unity 



at ~ ±45 deg (jWalkerl Il998l ). At lower latitudes still it will increase only slowly due to 
effects of refractive ISS. The increase in m pt between 90 and 45 deg partially compensates 
for the reduction in D(2d) due to increasing distance L. At still lower latitudes D(2d) might 
be expected to decrease. In the observations one sees a difference of low latitude behavior 
between the Northern and Southern hemispheres. This asymmetry can be understood by 
looking at the center panel of Figure [8] which shows that the emission measure is commonly 
higher for Southern latitudes. 

In a complete model the (unknown) compact fraction f c of the source flux density in 
the scintillating component must also be considered. Thus the expected variation of the 
ISS level with Galactic latitude must be combined with the probability distributions for the 
flux fraction and for the diameters of the compact components, which will further dilute the 
variation of D(2d) with latitude. 

In the Green Bank Interferometer observations cited above, the data were sampled daily 
(or on alternate days) over many years and so provided estimates of the actual modulation 
index, which were not reduced by the lengthening ISS time scale at low latitudes. However, 
even for these data the latitude dependence is not a strong effect. Note that the asymmetry 
about the Galactic plane is very similar to the asymmetry in the typical Ha emission as 
shown by the circles in the center plot. The lower panel of figure [8] plots the fraction of 
sources in the three timescale groups versus latitude and shows clear evidence that the fast 
scintillators dominate at high latitudes and that slow scintillators dominate at low latitudes. 
This agrees with the expected increase in ISS timescale at low latitudes outlined above. 



5.2. Dependence of ISS on Source Spectral Index and Flux Density 

As stated earlier the ISS of extragalactic sources is expected to be strongly suppressed, 
relative to that of pulsars, by the smoothing effect of their larger angular diameters. Con- 
sequently, we expect that the more compact sources should show higher levels of ISS. Here 
we examine the influence of mean flux density (S) and spectral index a (S oc u a ), since we 
expect synchrotron emitting sources to be more compact for larger a and lower S, due to 
the effects of synchrotron self-absorption and inverse Compton losses. 



In iHeeschenl ( 11984j )'s initial survey of short-term variability of a large sample of both 
steep-spectrum and flat-spectrum sources he found that the flat-spectrum sources varied, 
"flicker", but the steep-spectrum sources did not. This can be understood as the steep- 
spectrum sources are dominated by optically thin synchrotron emission with low brightness 



-19- 



temperatures, while the flat spectrum sources are dom inated by synchrotron self -absorbed 
components with very high brightness temperatures (jScheuer fc Williamsl Il968l ). making 
them compact enough to show ISS. 



The MASIV sources were selected to have flat spectral indices a > —0.3 (ILovell et al. 



20031 ) and so are predominantly quasars with compact cores. However, it is possible that 
there are also some very compact galaxies in the sample. Figure [9] shows the spectral 
index distributions separately for the sources with the visual classification as variable or 
non-var iable. The spectral indices shown are t hose used to form t h e sample: 1.4 GHz 



NVSS flCondon et al. Il998h to 8.5 GHz JVAS flPatnaik et al.l 11992 : Browne et al.l Il998 



Wilkinson et all Il998l ) or CLASS jMvers et all Il995h flux densities. This shows a slight 



increase in the fraction of sources that are variable with increasing a, in agreement with the 
expectation that the flatter (and inverted) spectrum sources are more compact. Though in 
Figure [10] the mean D(2d) shows no significant trend with spectral index, the bottom panel 
shows a slight increase in the fraction of variable sources for a > 0. This is mostly due to 
an increase in the fraction of slow variables which constitute the largest timescale group. It 
is worth pointing out that the surveys from which the flux densities were drawn to obtain 
spectral index were not coeval. It is likely then that any change in sample properties as a 
function of spectral index will be blurred as many of the sources vary intrinsically. 

Turning to the influence of mean flux density, we first discuss the visual variability 
classification and then plot D(2d) against flux density. The selection of sources for the 
MASIV survey divided them into a high mean flux density group (strong) and a low mean 
flux density group (weak S < 0.3 Jy). As reported in Paper 1 there was a greater fraction of 
variable sources in epoch 1 from the low flux density group than from the high flux density 
group. Combining all four epochs the numbers of weak sources that varied in 0-4 epochs is 
(94, 46, 36, 39, 33) and for strong sources the numbers are (62, 45, 49, 22, 23). Thus there 
are significantly more 3 & 4 time variables among the weak sources than among the strong 
ones, though this trend is not supported in the 2 or 1 time variables. 

Figure [UJ shows D(2d) versus mean flux density. The center panel shows a clear down- 
ward trend with increasing flux density in the lowest three bins, while in the fourth bin it 
increases but with fewer sources the mean has a large error. The interpretation is an in- 
creasing angular diameter of the compact source components with increasing total mean flux 
density. Furthermore, the bottom panel shows a decrease in the fraction of fast and medium 
scintillators with increasing mean flux density. These are exactly the trends expected if their 
effective angular diameters are constrained by a maximum brightness temperature due to 
self-absorption or inverse Compton losses (9 cx (S7Tb) ' 5 ). 
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5.3. Source Models 

The foregoing analysis establishes that about half of the 443 compact flat spectrum 
radio sources in the MASIV survey show ISS at an rms level above 1% over times of 2 d. 
We now consider a simple model for the compact source structure based on Appendix A. 
Equation flA5[) gives an approximate relation between D(2d) and parameters of the source 
(f c and 9 src ) and of the interstellar medium (L and V). We proceed by assuming a basic 
model for the latter parameters L = 500 pc and V = 50 km s -1 and finding constraints on 
the source. 

Figure [T21 shows (solid) contours of D(2d) versus the source parameters: compact com- 
ponent diameter and compact flux fraction - defined in the observer's frame. Also shown are 
(dashed) contours of Tb/Sj y . The majority of sources are in the range 0.0004 < D(2d) < 0.01 
which, for sources of 0.1 Jy mean flux density, maps to maximum brightness temperatures 
10 12 K to 10 14 K. These figures require substantial Doppler factors in the AGN jets comparable 
with those estimated from VLBI. However, we note that the plot of Figure [T21 only provides 
a guide since it is based on 500 pc as the distance to an interstellar scattering screen. 

The distribution of scattering electrons along the line of sight to each source is likely 
to be much more complex and can extend from tens of pc to a few kpc. Scattering at tens 
of pc has been shown to be important for the rare rapid scintillators (IHV) and scattering 
from more than 1 kpc is responsible for the slower timescale ISS associated from sight lines 
with large emission measures. Since the implied angular diameters scale roughly with the 
distance, the uncertainty in L corresponds to an extremely large range in implied brightness 
temperatures. We note, however, the fast ISS sources in our sample are mostly scintillating 
at levels of only 1-5%, unlike the large amplitude variations in the well studied IHV sources 
B0405-385, B1257-326 and J1819+3845. This suggests that the nearb y scattering r egions 



responsible are relatively rare covering only a small fraction of the sky. lLazio et al.l (120081 ) 
discuss the relative importance of sparsely distributed "clumps of scattering material" and a 
more uniformly distributed interstellar scattering plasma, suggesting that the former could 
be more important for ISS and the latter for angular broadening in the ISM. These authors 
find minimum diameters of ~ 1 — 2 mas at 1 GHz, which they suggest is caused by interstellar 
scattering, which predicts 40 — 80/xas when scaled to 5 GHz. It will be important to use 
the full MASIV data set to re-examine these questions, but since our emphasis here is on 
presenting the data, we postpone them to a later paper. 
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5.4. Combined effects of Intrinsic Variations and ISS 

Here we examine the competing contributions that scintillation and intrinsic variability 
would potentially make to the measured lightcurves. A source at an angular diameter dis- 
tance Da undergoing intrinsic variations on a timescale r has an implied maximum intrinsic 
angular size: 

= 17.3V ( T - //as, (4) 



lOOdaysy \lGpc 

where T> is the Doppler factor. Further, a variation in flux density AS implies an observed 
brightness temperature of 

ViOOmJyy VlOOdaysy ViGpcy V ; 

When mapped into the emission rest frame the brightness temperature is then reduced by a 
factor D 3 (l + z)~ 3 , under the hypothesis of intrinsic variation. 

As an example consider a source that undergoes 100 mJy intrinsic fluctuations in 100 days, 
as observed between epochs for several of our sources. At a typical distance Da ~ lGpc the 
implied maximum source size would be ~ 17T>/ms. Suppose further that it does not show 
ISS within an epoch, which implies it must be larger than 80 /xas in the observed frame and 
so T> > 4.6. Hence mapping equation (jSJ) into the emission frame gives TB )em ~5 x 10 10 K. We 
conclude that sources showing intra-epoch (intrinsic) variation and no ISS have relatively low 
emission frame brightness and, conversely, higher brightness sources that show intra-epoch 
variation have to show ISS. 



5.5. Dependence of ISS on Source Redshift 

We found redshifts for about half of the 443 sources i n the survey from th e published 



literature and we have subsequently measured another 69 (IPursimo et al.l 120081 ) for a total 



of 275 redshifts. This constitutes the largest sample of ISS measurements versus redshift. 

Figure [13] plots D(2d) versus redshift and reveals a highly significant decrease in the 
prevalence of ISS as redshift increases. In particular the middle panel shows that when 
binned in redshift the mean level of ISS drops steeply above redshift 2. As in other plots in 
this format the binned averages are independent of any SF threshold. However, note that 
the exact value of the lowest binned averages of D(2d) depend on the details of the noise 
subtraction. We note that, since the effect of confusion is to slightly increase D(2d), our 
estimates become upper limits when below the threshold of about 0.0004. 
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The bottom panel suggests that the fraction of fast variables drops more quickly with 
redshift than the fraction of slow and medium variables. However, the error bars show that 
this is only marginally significant, and in view of the importance of this question we list the 
source counts in Table If true it would imply that the drop in mean ISS level is due to an 
increase in angular diameter with redshift, which also lengthens the ISS timescale. While 
the drop in ISS level seen in the middle panel could be due either to an increase in diameter 
of the compact core of emission from these sources or a decrease in the compact fraction 
of flux density in that core, the latter interpretation would not explain a steeper decrease 
in fast ISS with redshift than in slow ISS. The most likely conclusion from this analysis is 
that the extremely compact emitting regions responsible for the ISS in over half the quasars 
studied appear broader in angular diameter with redshift above 2. The interpretation of 
this result involves either an evolution in the emitting regions with redshift or an angular 
broadening phenomenon due to propagation. We caution that a full consideration of selection 
effects must be made when interpreting this result. For example, redshifts are currently 
more complete in the strong (84%) than the weak (43%) sub-samples. Therefore, although 
this new cosmologically important result is the major finding from the M ASIV survey, we 



postp one a full discussion of the interpretation to a paper in preparation (jMacauart et al. 



20081 ) pending a thorough investigation of selection effects (IPursirno et al 



2008). Interested 



readers can co nsult preliminary d iscussions of the interpretation by lRickett et al.l (120071 ). We 
also note that lLazio et al.l (120081 ) plotted angular diameter at 1 GHz against redshift from 
a much smaller sample of scintillating and non-scintillating sources but could draw no firm 
conclusions. 



5.6. Intermittent Variability 

We now ask whether the sources that were only classified as variable in one to three 
of the four epochs are varying episodically or are the result of statistical uncertainty and a 
fixed threshold for variability on the raw modulation index. As a result of the measurement 
uncertainties there can be both false positives and false negatives, whose probabilities we 
can estimate. Concentrating on the 91 sources that were classified as variable in only one 
epoch and correcting for false positives leaves us with 61 sources ~ 13% of the total. These 
may be a category of short-lived, episodic scintillators revealed by our regular sampling over 
a full year. 



We note the strong case fo r intermittent ISS in the rapid variable PKS 0405-385 
( IKedziora-Chudczer et al.ll2001bl ). and so consider a simple model for the intermittency in 
terms of the longevity of each IDV episode and the duration between episodes. Intermittency 
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may arise either due to fluctuations in the level of turbulence responsible for ISS, or it may 
arise if the lifetimes of the bright, microarcsecond source components that undergo ISS are 
short. In the latter case one might expect sources to undergo IDV in conjunction with each 
outburst of the central engine. 

Consider a simple model in which the IDV episodes have a finite duration AT and an 
interval, Tf, between outbursts. Obviously in reality both quantities will have a distribution 
of possible values, but given the infrequency of our time sampling we restrict ourselves to this 
simple assumption here. For any given source the IDV episode commences at some random 
time U G [0, Tf, with the probability distribution of episodes distributed uniformly: 

p{U) = Tf\ 0<U< T f . (6) 

If we make a single observation the probability that the source will be exhibiting IDV is 

Pl = AT/7}, (7) 

so for a survey that examines N intermittent sources the expected number showing IDVs 
at any one time is NAT/Tf. One can similarly calculate the probability of detecting IDV 
in a source during one or more of multiple observing epochs of a multi-epoch survey (see 
Appendix [B]). In particular, the probability of detecting IDV in a source in one or more 
epochs of a four epoch survey, with epochs separated evenly in time by t Q b s , is 

{ AT+3f obs J- <-- /\rp 

> AT. (8) 

The number of IDV sources detected is a maximum when the interval between observing 
epochs exceeds the duration of IDV episodes because for shorter epoch intervals, after one 
merely discovers few new IDVs after the first epoch, one only keeps reobserving all the IDVs 
that were present in the first epoch. (Obviously, in the limit when t b s is small, multiple 
observations discover the same number of sources as a single-epoch survey.) 

Now, the mean detection rate of IDV sources in each epoch is 30%, whereas the fraction 
of sources that exhibited IDV in one or more of our four epochs is 58%. These two numbers 
imply a typical burst duration AT = 1.2 yr and a burst period of Tf = 3.8 yr. 

We can also calculate the corresponding probability that IDV is observed in a source in 
all four epochs: 

taH 0, ' 3tote >AT. 
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Based on the foregoing estimates of AT and Tf one estimates that only 4% of all sources 
should be common to all four epochs. However, the actual detection fraction is 12%. 

It should be remembered that this model does not take into account several effects 
which are likely to be important, including: (i) annual cycle effects influence the number of 
sources one detects at any one epoch (ii) there is likely a distribution of episode durations 
and repetition rates, (iii) the repetition is likely irregular and (iv) not all IDVs are likely to be 
episodic. We favor (ii), the importance of which is illustrated by the fact that many one-time 
IDVs were seen in epochs 2 & 3 that were (obviously) not detected in epoch 4. However, our 
model implies that sources that commenced IDV in these epochs should have been detected 
subsequently because the predicted burst duration exceeds the interval between observing 
runs (i.e. AT > 3t obs )- 

6. Conclusions 

We have reported results from the four epochs of the MASIV survey. There were 710 
sources with flat spectra (a < —0.3) near 5 GHz selected in weak and strong flux density 
groups surveyed for variability in four epochs over a year. These flat spectrum sources are 
predominantly quasars with compact emission probably from a core and jet, many with 
effective diameters small enough to show interstellar scintillation (ISS). In each epoch the 
flux density was measured using sub-arrays of the VLA every 2 hrs for about 12 hrs each 
day for 3-4 days. Sources were removed from the study if they showed evidence for changing 
correlated flux density due to confusion or resolution of their more extended structure, leaving 
443 sources which were analyzed for variability in two ways. 

The first was a binary classification based on the raw modulation index (visual method) 
in which 43% of the sources were classified as variable in 2, 3 or 4 of the epochs. The second 
was a fit to the epoch- averaged structure function parameterized by T>(2d) and a timescale 
r c h ar . In view of the uncertainties in the latter we classified sources as slow (> 3 days), 
medium (0.5-3 days) or fast (< 0.5 days) if T>(2d) exceeded 4 x 10~ 4 . By this criterion 37% 
of the sources varied with more than 1.4% modulation index over 2 days, which is similar to 
the 43% variables by the visual classification. 

We found that _D(2d) and timescale varied both with coordinates in the Galaxy and also 
with source-based quantities. This confirms that the variations are dominated by ISS, which 
depends on both the strength of scattering and the distance to the scattering region and 
also on the fraction of flux density in its most compact component and its effective angular 
diameter. The following is a summary of our findings: 
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The amplitude of 2-d variability increases with increasing emission measure estimated 
from Ha intensity for each line of sight. Emission measure is the column density for 
the square of the electron density which is expected to be strongly correlated with in- 
homogeneity in the ionized medium that causes ISS. This result provides observational 
evidence that ISS is the dominant cause of the variations. We find fast variations dom- 
inate for low emission measure, as expected since such regions will be seen out of the 
plane and closer to the Earth, and slow variations dominate for high emission measure 
which are typically seen at greater distances toward the Galactic plane especially for 
southern latitudes where the Ha intensity is high. 

The amplitude of 2-d ISS variability varies significantly with Galactic latitude but 
differs substantially between positive and negative latitudes. The expected behaviour 
is complicated; greater path lengths at low latitudes, where the scattering should be 
stronger, cause the scattering to be slower which should reduce the rms over 3 d. 
However, the observed timescales show that there are more sources with fast varia- 
tions at high latitudes and more sources with slow variations at low latitudes in both 
hemispheres, in clear support of ISS as the dominant cause. 

The ISS modulation index tends to decrease with increasing mean flux density, as 
expected if the compact emission is limited by synchrotron self absorption or inverse 
Compton losses to have a maximum brightness temperature. In that case the expected 
angular diameter oc \/S which will quench the ISS of the stronger sources. 

There is little change in the ISS amplitude with spectral index for our sample with 
a > -0.3. 

There is evidence that the ISS can be intermittent on times of 3-6 months for some 
sources, but this is hard to quantify from the 3 day observing sequences, when the 
time scale of the variations is of the same order. 

We model D(2d) as a function of compact source component fractional flux density 
and angular diameter, from which we find compact diameter to lie in the range 0.005 
- 0.15 milli arcseconds and brightness temperatures in the range 10 12 — 10 14 K. 

The most far-reaching result reported here is the discovery of a decrease in both the 
fraction of sources that scintillate and in their scintillation amplitude beyond redshifts 
around 2. We conclude that there is an increase in the typical angular diameter of 
the most compact radio-emitting regions of the quasars beyond reshift 2. The pos- 
sible interpretations o f this exciting result will be presented in a companion paper 



flMacquart et all 12008 ). 
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A further surprise (at least to us) was the apparent absence of the very rapid variables 
(IHV). J1819+3845 fell in our sample, but it was the only source to show remarkable 
large a mplitude variabili ty. J0929+5013 showed rapid variability in the January 2002 
epoch (ILovell et al.ll2003l ) but, although monitored closely, revealed only slower, many- 
hour variability in the three later epochs. We had expected to find more of these rapid 
variables especially given that two of the three known, J1819+3845 and PKS 1257-326 
were discovered serendipitously. This strongly suggests that the IHV sources lie behind 
discrete local interstellar clouds which cover a small fraction of the sky. 
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A. Structure function for ISS 

In all of the MASIV data the intrinsic source diameters (8 src ) are large enough to suppress 
the scintillations relative to those of a point source (such as a pulsar), which at 6 cm would 
be expected to have a true modulation index m pt near unity. In this section we describe a 
model for the structure function for such extended sources; this is easier to interpret than 
the apparent modulation index mj. 

The source diameter smoothes the ISS of a point source and so reduces the modulation 
index to tuq and lengthens the timescale to tq. If the scattering is concentrated at a distance 
L from the Earth and we are near the transition from weak to strong scintillation, a useful 
approximate relation is: 







7/6 



m 



P t 



(^ + C) 7/12 ' 



(Al) 



Here 9-p = a/ A/ (2tyL) is the angular size subtended by the Fresnel scale (rp) and 6 src is the 
angular radius of the sour ce. The exponents 7/6 and 7/12 apply for a Kolmogorov spectrum 
of interstellar turbulence (jColes et al.lll987l ). To the same accuracy the ISS timescale for a 
point source would be rp = LQ-p/V and the source smoothing would increase it to: 



Tg = Tp- 



(A2) 



where V is the effective transverse velocity of the observer relative to the layer of scattering 
plasma. Note that when the source diameter is sufficiently large to suppress the scintillations 
we have 



m e ~ (Of/8, 



, 7/6 Ta 
src I i'd 



L8 SIC /V . 



(A3) 



where we have set m pt = 1 (jRickettlll986l ) 



The theoretical form of the auto-correlation function for an extended source that sub- 
st antially suppresses the scintillation index is given by the low wavenumber approximation 
of IColes et al.l (119871 ). This in turn gives the theoretical structure function, whose detailed 
shape depends on both the spectrum of the plas ma density and on it s distribution along 
the line of sight through the Galaxy. Figure 14 of iRickett et al.l (120061 ) shows the form for 
sources with Gaus sian brightness with pea k temperature 10 11 — 10 13 K, when the medium 
is modeled by the ICordes and Laziol (120051 ) electron distribution at a Galactic latitude of 
45° (away from the Galactic center). The form of the structure function at small time lags 
depends strongly on the density distribution in the local ISM. A useful approximation to 
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these results is given by: 



m = Vlrnl-^— a , (A4) 



where 1 < a < 2 is a constant that depends on the density distribution in the local ISM. Here 
a ~ 2 for a local bubble with low turbulence such that the effective scattering distance is 
beyond the bubble (> 100 pc) and alternatively a ~ 1 if the medium is uniformly turbulent 
out to a scale height (as in the "disk" of the CL05 electron density model). We have also 
introduced an extra variable f c that is the fraction of the source flux density in the bright 
(core) component. 

Equations (1A4I) and (1A3I) thus provide a simple interpretation for our estimates of D{t = 
2d). Inserting me and t$ from equation (IA3I) we obtain: 

An important feature of this result is that D(2d) decreases steeply with increasing 8 src and so 
provides a sensitive measure of source diameter. In estimating D(2d) for every source-epoch 
we set a = 1, which allows us to estimate the scintillation timescale tq without having to also 
estimate the exponent a. D(2d) can be converted to an effective 2-day modulation index 
by m 2 d = A/0.5-D(2d). We note that m 2( j can substantially exceed the apparent modulation 
index mj when the time scale is longer than about 2 days. 



B. A simple model for IDV intermittency 

Consider a model in which a source outbursts every duration Tf in time and each 
IDV episode lasts AT. The initial outburst time is unknown, but its probability is evenly 
distributed in the interval ^ G [0,X/): p{ti) = TJ . Now consider a function ]{t) = 1 — 
[H(t) — H(t — AT)], which assumes the value one whenever the source shows no IDV but 
the value zero when it is on. 

Thus the fraction of the time in which the source is off for initial burst durations between 
ti and ti + d8ti is f(ti)p(ti)5ti. Thus the probability that IDV is absent is 

fTf _ 

Piofr = / dtif(ti)p(ti) 







AT , . 

= 1-^- (Bl) 

Thus the probability that the source is observed to exhibit IDV is 1 — p\ g = AT/Tj. 
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Now suppose we look for IDV at times t — 0, t bs, 2t bs, 3t bs- The fraction of the burst 
times between ti and ti + 5ti for which IDV is absent in all four observations is f{ti)f(ti + 
t bs) f(ti + 2t ohs ) f(ti + 3t b s )p(ti)5ti. If we assume that the repetition period exceeds the total 
duration of our observations (i.e. Tf > 3t b s ), the probability of observing no IDV over all 
four epochs is 

r T s 

P4off = / dt i f(t i )f(t i + t ohs )f(t i + 2t obs )f(t i + 3t ohs )p(t i ) 
Jo 

1 _ Ar+3t obs . ^ \rp 

AAT f ' ' ( B2 ) 

1 jT~, tobs > AT. 

Thus the probability that IDV is observed in any one or more of these four epochs is 

- i _ - / ^f^' tohs ~ AT ' 

P&ny of 4 on _ 1 Pioff — ^ 4AT t h > AT 

We can similarly consider the probability of observing IDV in multiple observations by 
defining the function f(t) = H(t) — H(t — AT) which takes the value one whenever the IDV 
is on and zero otherwise. The probability that IDV is observed in all four epochs is thus 



£>4on = 



dtif(ti)f(ti + t ohs )f{U + 2t obB )/(t i + 3t obB )p(ti) 

AT-3t obs OJ- <- \rp 

0, T f > 3t ohs > AT. 
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Epoch 1 Epoch 1 




<S> (Jy) <S> (Jy) 



Fig. 1. — Left: Raw modulation index for epoch 1 plotted against mean source flux density. 
Sources classified as variable plotted as circles and non-variable as pluses. Right: Similar 
plot for rri2h as described and defined in §4.11 In both plots the dashed line is equation(OQ) 
with s = 0.0015 Jy and p = 0.01 and the solid line has s = 0.0013 Jy and p = 0.005. Similar 
plots are obtained for the other three epochs. 
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Fig. 2. — Upper: Time series of flux density (Jy) for J2325+3957 versus day number from 
2002 January 1; error bars are as in Paper 1. Lower: Structure Function D(t) of flux density 
(normalized by its mean) averaged over all 4 epochs. The dashed line is the estimated noise 
level -D n oise and the solid line is a simple model fit (see §4.31) . The vertical bar centered near 
lag of 2 d is an estimate of .0(2 day) and its standard error. 
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Fig. 3.— Examples of MASIV variability. J0949+5819 (left) and J1328+6221 (right). The 
top panel in each case shows all the flux densities in Jy against day number from 2002 
January 1. The middle four panels show light-curves for each of the first four epochs. The 
horizontal scale is the same (four days) in each case. Lower panels show structure function 
of flux density (normalized by its mean) averaged over all 4 epochs. Dashed line is the 
estimated noise level Aioise and the solid line is a simple model fit (see text). J0949+5819 
shows evidence of episodic scintillation, possibly related to intrinsic source changes as there 
appears to be a corresponsence between mean flux density over an epoch and amplitude 
of scintillation. J1328+6221 shows strong, consistent scintillation over all four observation 
epochs. 
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Epoch 1 




<S> (Jy) 

Fig. 4. — = y/D(2 d) (without any noise subtraction) plotted against mean source flux 
density. Circles are sources classified as variables and pluses as non-variables. The lines are 
equation^ with s = 0.0015 Jy and p = 0.01 for the dashed line and s = 0.0013 Jy and 
p = 0.005 for the solid line. Similar plots are obtained for the other three epochs. 
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1 2 3 >3 

Number of Inflection Points 




2002 Jan 2002 May 2002 Sep 2003 Jan 



Epoch 

Fig. 5. — Upper: The number of sources classified as variable versus the observed number of 
changes in the sign of the derivative of flux density vs time (i.e. number of inflection points). 
Clearly a majority of the sources vary on timescales of 3 days or longer. Lower: Numbers of 
sources classified as variable in each epoch. Screens moving at the LSR would be expected 
to result in fewer variables beeing seen in epoch 3; we see no such deficit. 
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Fig. 6. — The number of "variable" epochs plotted against the value of _D(2d) fitted to the 
cumulative structure function for all four epochs. The large circles show the mean values for 
each group of sources (see Table [3]). 
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10 10 
WHAM Ha Intensity (Rayleighs) 



Fig. 7. — Upper panel: Log-log Scatter plot of D(2d) against WHAM Ha emission which 
is proportional to the emission measu re on a line of sight sampled on a 1 degree grid of 
the Northern sky (IHaffner et al.ll2003l ). The intra-hour variable source J1819+3845 is off 
scale at 0.25. The different symbols represent the three classifications of ISS timescale, as 
described in the text. Center panel: Mean value of D(2d) in the indicated bins of Ho; 
emission including all sources except J1819+3845; vertical bar gives the standard error in 
the mean. Lower panel: Fraction of sources above the threshold in each timescale class in 
each bin showing that fast ISS is commonest for the lower column density of electrons and 
slower ISS dominates for higher column densities. Error bars assume binomial distributions. 
The same method is used in Figures [TUl [IT] & [131 
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Fig. 8. — Upper panel: Scatter plot of _D(2d) (on log scale) against Galactic latitude. The 
intra-hour variable source J1819+3845 is off scale at 0.25. The different symbols represent 
the three classifications of timescale, as described in §4.31 Center panel: Squares show 
mean value of _D(2d) in 30 degree bins of latitude including all sources except J1819+3845; 
vertical bar gives the error in the mean. Note the N-S asymmetry in the circles show mean 
Ha emission in each bin. Lower panel: Fraction of sources above the threshold in each 
timescale class in each bin. 
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Fig. 10. — D(2d) as a function of spectral index. Upper panel: Logarithmic scatter plot of 
D(2d). Different symbols represent the three classifications of ISS timescale. Center panel: 
Mean D(2d) binned by spectral index including all sources except J1819+3845. Lower panel: 
Fraction of sources above the threshold in each timescale class in each bin. No significant 
trend is seen. 
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Fig. 11. — Upper panel: Scatter plot of D(2d) against mean flux density. The higher values 
are more common for lower flux density sources. The different symbols represent the three 
classifications of ISS timescale, as described in the text. Center panel: Mean value of D(2d) 
in bins of mean flux density for all sources (excluding extreme IHV quasar J1819+3845). 
Note lower levels of ISS for the sources with higher flux density. The points centered at 2.5 
Jy have only a few sources in each timescale group giving larger errors in the mean for that 
bin. Lower panel: Fraction of sources in each timescale class in each bin. Note decreasing 
occurrence of fast variables among stronger sources. 
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Fig. 12. — Solid lines show contours of log 10 [D(2d)] versus f c and 9 src based on a Kolmogorov 
model for ISS in a region at a distance of 500 pc. Dashed lines are contours of log 10 [T b /Sj y ] 
for total source flux density Sj y . 
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Fig. 13. — Upper panel: Scatter plot of _D(2d) against source redshift. The different symbols 
represent the three classifications of ISS timescale, as described in the text. Center panel: 
Mean value of _D(2d) in redshift bins for the 271 sources (out of 443) with measured redshift 
(excluding extreme IHV quasar J1819+3845). Note lower levels of ISS at high redshift. 
Values below the dash-dot line are upper bounds since they may be raised slightly by low 
level confusion. Lower panel: Fraction of sources in each timescale class in each bin. The 
two bins above z=3 have been combined, as there are so few sources at this redshift. 
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Table 1: A Summary of the observation dates, durations and VLA array configurations for 
the first four MASIV epochs. 



Epoch Number 


Start Time (UT) 


Duration (h) 


Array Configuration 


1 


2002 January 19,19:14 


72 


D -> A 


2 


2002 May 9,08:02 


72 


A -> BnA 


3 


2002 September 13, 17:50 


96 


B -> CnB 


4 


2003 January 10, 16:21 


72 


C -> DnC 



Table 2: Observed numbers and percentages of sources classified as variable in the four 
epochs. Percentages are given of the 443 sources observed, which excludes those used as 
calibrators in more than one epoch. The predicted percentage of misclassifications assuming 
all sources were non- variable (where P = 0.05 is the probability of a single misclassification). 



Number of 
variable epochs 


Observed number 
of sources 


Observed percentage 
of sources 


Predicted fraction 
of misclassifications 


Predicted fraction 
if none are variable 





161 


36% 


(1-P) 4 


81.5% 


1 


90 


21% 


AP(l-Pf 


17.2% 


2 


79 


18% 


6P 2 {1-P) 2 


1.3% 


3 


58 


13% 


4P 3 (1 - P) 


0.05% 


4 


55 


12% 


p4 


0.0006% 
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Table 3: Statistics of SF and visual variability classifications for the 443 source sample. 



No. epochs "variable" 


No. of Sources 


% Sources 


Mean D(2d) 


No. > 4 x 10~ 4 


No. < 4 x lO" 4 





161 


36% 


0.00024 


27 


134 


1 


90 


21% 


0.00071 


54 


36 


2 


79 


18% 


0.0011 


64 


15 


3 


58 


13% 


0.0027 


58 





4 


55 


12% 


0.0071 


55 






Table 4. MASIV source variability characteristics. Column 4 gives Ha in Rayleighs for WHAM. Columns 5-8 give 
mean flux density for epochs 1-4; Columns 9-12 give raw modulation index in % for epochs 1-4; Column 13 a err as %; 
Column 14 gives the visual classification per epoch for variable (y), non- variable (n) and non- variable source used as a 

secondary calibrator (c); Column 15 gives D(2d) from fit to epoch average SF where error -1 flags a non- convergent 
fit; Column 16 gives timescale classifications, l:r cnar < 0.5d, 2:0. 5d< r char < 3d, 3:r char > 3d, 0: -D(2d) below threshold 

of 4 x 1CT 4 . 



Name RA(J2000) Dcc(J2000) Ha R S„(Jy) Mod index % %err Var D(2d) tim 



J0005+3820 


00 


:05 


:57 


.17 


38 


:20 


:15. 


1 


2.9 


0. 


61 


0. 


.57 





.43 


0. 


.47 


1 


.1 


1 


.1 


1 


.2 


1 


.5 





.57 


nnny 


+4. 


73e-04 


± 


1.4e-05 


3 


J0007+5706 


00 


:07 


:48 


.47 


57 


:06 


:10. 


4 


10.8 


0. 


18 





20 





.17 





.17 


2. 


.1 


0. 


9 


1 


.1 


1 


.0 





.88 


nnnn 


+2. 


91c-04 


± 


4.3c-05 





J0009+1513 


00 


:09 


:43 


.47 


15 


:13 


:37. 


9 


0.7 


0. 


16 


0. 


16 





.17 





16 


1 


.0 


1 


.7 


2 


.0 


1 


.3 





.94 


nnnn 


+4. 


92e-04 


± 


5.5e-05 


2 


J0010+1724 


00 


:10 


:33 


.99 


17 


:24 


:18. 


.7 


0.8 


0. 


85 





83 





.81 


0. 


81 





.8 


1 


.0 





.9 





.5 





.53 


nnnn 


+1. 


07c-04 


± 


1.3c-05 





J0010+4412 


00 


:10 


:30 


.04 


44 


:12: 


:42. 


.4 


2.3 


0. 


20 





20 





.19 





.18 


1. 


.5 


1 


,2 


1 


.5 


1 


.7 





.83 


nynn 


+4. 


28e-04 


± 


3.9e-05 


1 


J0011+7045 


00 


:11 


:31 


.90 


70 


:45 


:31. 


6 


6.7 


0. 


62 





.71 





.63 





.59 


1 


.9 





6 





.6 


1 


.4 





.54 


ynny 


+3. 


46e-04 


± 


3.6e-05 





J0017+5312 


00 


:17 


:51 


.75 


5.3 


:12 


:19. 


1 


13.3 


0. 


61 





60 





.61 


0. 


.65 


1 


.5 





.9 


1 


.2 


1 


.2 





.54 


nnyy 


+3. 


88e-04 


± 


-l.Oc+00 





J0017+8135 


00 


:17: 


:08 


.47 


81 


:35: 


:08. 


1 


2.2 


0. 


91 





92 





.92 





.92 


1 


.4 


1 


.2 





.5 





.7 





.52 


nnnn 


+2. 


13e-04 


± 


1.6e-05 





J0019+2021 


00 


:19 


:37 


.85 


20 


:21 


:45. 


5 


0.7 


1. 


16 


1 


05 


1 


.03 


1 


22 


2 


.0 





.2 


1 


.0 





.7 





.52 


ynnn 


+4. 


86e-04 


± 


3.6e-05 


3 


J0019+7327 


00 


:19 


:45 


.78 


73 


:27 


:30. 





4.3 


0. 


45 





46 





.47 





.47 


1 


.5 





.8 


1 


.7 


1 


.0 





.57 


nnyn 


+5 


27c-04 


± 


1.0c+03 


3 


J0038+4137 


00 


:38 


:24 


.84 


41 


:37 


:06. 





2.6 





53 





53 





.53 





53 


1. 


.3 





.9 


1 


.6 


1 


.0 





,56 


ynyn 


+5. 


71c-04 


± 


2.6e-05 


3 


J0041+1339 


00 


:41 


:17 


.21 


13 


:39 


;27. 


4 


1.1 


0. 


15 





16 





.16 





16 


0. 


.8 


2 


.5 


1 


.9 


1 


.8 





.97 


nynn 


+7. 


29e-04 


± 


5.7e-05 


2 


J0042+2320 


00 


:42 


:04 


.54 


23 


:20 


:01. 


1 


0.7 


0. 


.78 





79 





.84 


0. 


.87 





.7 


1 


.0 





.8 





.1 





.52 


nnnn 


+1. 


87c-04 


± 


1.0c+03 





J0042+5708 


00 


:42 


:19 


.45 


57 


:08 


:36. 


5 


23.5 


0. 


86 


0. 


84 





.86 


0. 


.83 


1 


.4 





.9 





.4 





.5 





.52 


ynnn 


+1. 


06c-04 


± 


9.5e-06 





J0043+0502 


00 


:43 


:46 


.73 


05 


:02 


:56. 





0.4 


0. 


.10 





10 





.10 





10 


1 


1 


2 


.3 


2 


.2 


1. 


.8 


1 


.40 


nnnn 


+7. 


44e-04 


± 


8.5e-05 


3 


J0052+1633 


00 


:52 


:36 


.16 


16 


:33 


:00. 


4 


0.8 


0. 


17 





17 





.17 





.17 





.8 


1 


.0 


1 


.1 





.4 





.91 


nnnn 


+1. 


39e-05 


± 


1.0c+03 





J0056+1625 


00 


:56 


:55 


.29 


16 


:25 


:13. 


3 


0.8 


0. 


19 


0. 


23 





.29 


0. 


.27 


1 


.7 


3 


.0 


2 


.9 


1 


.8 





.72 


nyyy 


+1. 


39e-03 


± 


l.le-04 


1 


J0056+3142 


00 


:56 


:22 


.61 


31 


:42 


:09. 


1 


1.2 


0. 


08 





.08 





.07 





07 


2 


.3 


2. 


.2 


2 


.2 


2 


.6 


1 


.82 


nnnn 


+5. 


83e-04 


± 


9.3e-05 


1 


J0057+3021 


00 


:57 


:48 


.88 


30 


:21 


:08. 


.8 


1.2 


0. 


67 





69 





.70 


0. 


.67 





6 


1 


.0 


1 


.3 


1 


.2 





.53 


nnnn 


+2. 


36e-04 


± 


2.4e-05 





J0101+5004 


01 


:01 


:16 


.99 


50 


:04 


:44. 


9 


7.0 


0. 


16 





16 





.17 


0. 


16 


2 


.6 


1 


.1 


1 


.7 


2 


.9 





.95 


ynny 


+1. 


16e-03 


± 


l.le-04 


2 


J0102+5824 


01 


:02 


:45 


.76 


58: 


:24 


:11. 


1 


17.7 


1. 


37 


1 


88 


2 


.63 


2. 


34 


8. 


.0 


3 


.5 


6 


.3 


3 


.1 





.53 


yyyy 


+1. 


09e-02 


± 


3.4e-04 


3 


J0104+1134 


01 


:04 


:10 


.62 


11 


:34: 


:16. 


3 


0.8 


0. 


.08 


0. 


.08 





.08 





08 


1 


.4 


2 


.5 


1 


.8 


1 


.7 


1 


.70 


nnnn 


+3. 


64e-04 


± 


8.5e-05 





J0104+2949 


01 


:04 


:08 


.32 


29 


:49 


:05. 


6 


1.2 


0. 


09 





09 





.09 





.09 


2 


.8 


3 


.5 


2 


.2 


1 


.5 


1 


.54 


yynn 


+9. 


06e-04 


± 


1.4e-04 


1 


J0105+4819 


01 


:05 


:49 


.92 


48 


:19 


:03. 


1 


5.4 


1. 


09 


1 


17 


1 


.15 


1 


.19 





.8 





.2 


1 


.0 





.2 





.51 


nnnn 


+1. 


87c-04 


± 


2.2e-05 





J0108+0135 


01 


:08 


:38 


.77 


01 


:35 


:00. 


3 


0.7 


2. 


.41 


2 


.41 


2 


.60 


2. 


.83 





.5 


1 


.7 


1 


.3 





.8 





.50 


nyyn 


+2. 


56e-04 


± 


3.6e-05 





.10111+3906 


01 


:11 


:37 


.31 


39 


:06 


:28. 





1.8 


1. 


31 


1 


.34 


1 


.29 


1. 


.31 





.2 





.6 





.7 





.7 





.51 


nnnn 


+4. 


77c-05 


± 


6.6e-06 





.10121+0422 


01 


:21 


:56 


.86 


04 


:22 


:24. 


7 


0.6 


0. 


89 





.89 





.95 





.97 


1 


6 


2 


.0 


1 


.2 





.7 





.52 


yynn 


+5. 


13e-04 


± 


2.9e-05 


2 


J0121+1149 


01 


:21 


:41 


.59 


11 


:49 


:50. 


4 


0.7 


1. 


70 


1 


.71 


1 


.57 


1. 


59 





.7 


1 


.4 





.5 





.8 





.51 


nnnn 


+1. 


34e-04 


± 


1.8e-05 





J0122+0310 


01 


:22 


:01 


.91 


03 


:10 


:02. 


4 


0.5 


0. 


12 





.11 





.10 





.10 





.7 


1 


.6 


1 


.6 


1 


.4 


1 


.30 


nnnn 


+1. 


30e-04 


± 


5.1e-05 





J0122+2502 


01 


:22 


:38 


.81 


25 


:02 


:31. 


.7 


0.9 


0. 


82 





.82 





.84 


0. 


.86 


1 


1 


1 


.5 


1 


.1 


1 


.2 





.52 


nnny 


+3. 


81c-04 


± 


2.5e-05 





J0122+2536 


01 


:22 


:40 


.77 


25 


:36: 


:14. 


9 


0.9 


0. 


10 


0. 


11 





.11 





12 


1. 


.5 


2 


.0 


1 


.8 


4 


.8 


1 


.27 


nnny 


+2. 


16e-03 


± 


3.9e-04 


1 


J0126+2559 


01 


:26 


:42 


.79 


25 


:59 


:01. 


2 


1.0 


0. 


.74 





.76 





.79 


0. 


.81 





.8 


2. 


.1 





.8 





.5 





.53 


nynn 


+3. 


68e-04 


± 


6.2e-05 





J0135+2158 


01 


:35 


:53 


.49 


21 


:58 


:16. 


1 


0.9 


0. 


13 


0. 


.13 





.12 


0. 


12 


1 


.5 


1 


.6 


1 


.5 


1 


.1 


1 


.15 


nnnn 


+3. 


04e-04 


± 


2.7c-05 





J0136+4751 


01 


:36: 


:58 


.59 


47 


:51 


:29. 


1 


8.0 


2. 


75 


3 


26 


2 


.52 


2. 


.71 


1 


.8 


2 


.0 


1 


.2 


4 


.2 





.51 


yyyy 


+2. 


23e-03 


± 


6.4e-05 


3 


J0141+7506 


01 


:41: 


:00 


.62 


75 


:06 


:28. 


2 


1.9 





09 


0. 


10 





.09 





09 


4. 


.3 


1 


.4 


2 


.1 


2 


.0 


1 


.49 


ynnn 


+1. 


72e-03 


± 


-1.0e+00 


3 



Table 4 — Continued 



Name 


RA(J2000) 


Dec(J2000) 


Ha R 




S„(Jy) 




Mod index % 


%err 


Var 


D(2d) 




tim 


J0150+2646 


01:50:02.80 


26:46:28.0 


1.3 


0.12 


0.14 0.13 


0.12 


7.7 


4.2 


4.2 


7.2 


1.14 


yyyy 


+9.10c-03 


± 


1. 


9c-03 


1 


J0151+2744 


01:51:27.14 


27:44:41.7 


1.4 


0.85 


0.84 0.83 


0.83 


0.2 


0.8 


0.8 


0.4 


0.52 


nnnn 


+9.88e-05 


± 


2. 


Oc-05 





J0152+2207 


01:52:18.05 


22:07:07.7 


0.9 


1.06 


1.04 1.04 


1.06 


1.2 


1.0 


0.7 


0.6 


0.52 


ynnn 


+1.25e-04 


± 


1 


5c-05 





J0154+4743 


01:54:56.29 


47:43:26.5 


8.4 


0.44 


0.35 0.39 


0.39 


1.3 


1.0 


2.0 


2.0 


0.60 


ynyy 


+8.34c-04 


± 


4. 


9c-05 


3 


J0158+4137 


01:58:23.84 


41:37:39.5 


2.6 


0.08 


0.08 0.08 


0.09 


1.9 


1.6 


1.4 


2.0 


1.63 


nnnn 


+l.llc-04 


± 


3 


4c-05 





J0200+0322 


02:00:40.81 


03:22:49.5 


0.5 


0.32 


0.27 0.22 


0.20 


1.0 


2.8 


1.2 


0.9 


0.75 


nynn 


+6.87c-04 


± 


8, 


3e-05 


1 


J0204+1514 


02:04:50.41 


15:14:11.0 


0.7 


2.81 


3.07 3.04 


3.01 


0.9 


1.3 


0.1 


0.8 


0.50 


nnnn 


+2.38e-04 


± 


2. 


9c-05 





J0212+3626 


02:12:09.73 


36:26:17.9 


2.2 


0.09 


0.08 0.08 


0.08 


2.3 


1.8 


1.4 


2.1 


1.68 


nnnn 


+3.78c-04 


± 


9 


3c-05 





J0213+3652 


02:13:48.19 


36:52:33.9 


2.1 


0.13 


0.16 0.16 


0.18 


1.6 


1.2 


2.0 


1.0 


0.96 


nnyn 


+5.35e-04 


± 


1 


Oc+03 


1 


J0217+7349 


02:17:30.81 


73:49:32.6 


2.2 


3.43 


3.46 3.41 


3.41 


0.1 


1.0 


0.1 


0.1 


0.50 


nnnn 


+3.95e-04 


± 


2. 


7c-05 





J0228+6721 


02:28:50.05 


67:21:03.0 


4.3 


2.07 


2.02 1.79 


1.63 


0.6 


0.1 


0.6 


0.3 


0.51 


nnnn 


+2.30c-05 


± 


4. 


7c-06 





J0237+2046 


02:37:50.62 


20:46:18.4 


1.4 


0.10 


0.10 0.10 


0.10 


3.5 


5.2 


2.8 


2.4 


1.41 


yyyy 


+4.10e-03 


± 


2. 


3c-04 


3 


J0237+2848 


02:37:52.40 


28:48:08.9 


1.8 


2.78 


3.05 3.47 


3.69 


0.9 


0.4 


1.1 


0.1 


0.50 


nnnn 


+3.08e-04 


± 


5 


5c-05 





J0238+1636 


02:38:38.93 


16:36:59.2 


1.1 


1.50 


1.61 1.58 


1.34 


2.8 


3.0 


1.2 


1.8 


0.51 


yyyy 


+1.19e-03 


± 


9 


5c-05 


1 


J0239+0416 


02:39:51.26 


04:16:21.4 


0.8 


0.65 


0.60 0.60 


0.68 


0.7 


0.9 


1.6 


0.5 


0.54 


nnnn 


+3.06e-04 


± 


4. 


.5c-05 





J0242+1101 


02:42:29.17 


11:01:00.7 


1.0 


1.12 


1.16 1.24 


1.21 


0.2 


1.1 


0.9 


0.8 


0.51 


nnnn 


+1.41c-04 


± 


2. 


2c-05 





J0246+3536 


02:46:21.07 


35:36:37.9 


3.0 


0.17 


0.17 0.16 


0.16 


1.6 


4.1 


2.1 


1.1 


0.93 


nyyn 


+2.20c-03 


± 


1. 


.5e-04 


2 


J0248+0434 


02:48:14.82 


04:34:40.8 


0.8 


0.15 


0.15 0.15 


0.14 


0.9 


1.0 


1.1 


1.3 


1.01 


nnnn 


+4.21c-05 


± 


1 


8c-05 





J0249+0619 


02:49:18.01 


06:19:51.9 


0.8 


0.68 


0.68 0.70 


0.68 


0.9 


1.1 


1.1 


0.9 


0.53 


nnnn 


+1.75e-04 


± 


2. 


3c-05 





J0251+4315 


02:51:34.53 


43:15:15.8 


5.7 


1.18 


1.23 1.20 


1.21 


0.6 


0.6 


0.2 


0.2 


0.51 


nnnn 


+6.50e-06 


± 


4. 


le-06 





J0253+3217 


02:53:33.65 


32:17:20.8 


2.7 


0.09 


0.11 0.12 


0.14 


2.5 12.5 6.6 4.3 


1.21 


yyyy 


+2.23e-02 


± 


1. 


5c-03 


3 


J0257+4338 


02:57:59.07 


43:38:37.6 


5.6 


0.13 


0.14 0.14 


0.14 


1.5 


0.8 


1.4 


1.2 


1.07 


nnnn 


+2.11c-04 


± 


2. 


4c-05 





J0259+1925 


02:59:29.65 


19:25:44.3 


1.3 


0.20 


0.23 0.22 


0.19 


1.6 


4.8 


2.4 


3.1 


0.80 


nyyy 


+2.65e-03 


± 


2. 


6e-04 


1 


J0301+0221 


03:01:00.08 


02:21:46.4 


2.3 


0.14 


0.13 0.14 


0.13 


1.2 


1.4 


1.4 


1.0 


1.09 


nnnn 


+8.83c-05 


± 


2. 


lc-05 





J0303+4716 


03:03:35.24 


47:16:16.2 


5.0 


1.28 


1.29 1.41 


1.22 


0.4 


1.1 


0.5 


0.4 


0.51 


nnnn 


+1.12e-04 


± 


9 


2c-06 





J0304+6821 


03:04:22.00 


68:21:37.4 


2.6 


0.65 


0.67 0.68 


0.67 


1.8 


1.4 


0.7 


0.7 


0.54 


yynn 


+3.29e-04 


± 


3 


7c-05 





J0308+1208 


03:08:39.38 


12:08:49.7 


1.9 


0.08 


0.08 0.08 


0.07 


4.5 


3.7 


4.0 


1.2 


1.76 


yyyn 


+3.26e-03 


± 


2. 


le-04 


2 


J0309+1029 


03:09:03.62 


10:29:16.3 


1.5 


0.65 


0.69 0.73 


0.71 


0.7 


1.2 


1.2 


0.5 


0.53 


nnnn 


+3.51c-04 


± 


2. 


2e-05 





J0313+0228 


03:13:13.40 


02:28:35.2 


6.8 


0.13 


0.13 0.15 


0.13 


6.3 


6.9 


6.7 


3.0 


1.07 


yyyy 


+1.19e-02 


± 


5. 


9e-04 


3 


J0314+4345 


03:14:08.05 


43:45:19.7 


3.5 


0.20 


0.21 0.20 


0.20 


1.1 


1.1 


1.1 


0.7 


0.81 


nnnn 


+1.35c-04 


± 


2. 


5e-05 





J0319+4130 


03:19:48.16 


41:30:42.1 


4.6 


18.10 18.10 17.40 17.30 


0.5 


0.5 


0.7 


0.1 


0.50 


nnnn 


+3.82c-05 


± 


7 


9c-06 





J0321+1221 


03:21:53.10 


12:21:13.9 


2.1 


1.50 1.52 1.51 1.49 


0.7 


1.6 


1.0 


0.8 


0.51 


nnnn 


+2.36e-04 


± 


3 


Oc-05 





J0323+0446 


03:23:14.72 


04:46:12.5 


4.6 


0.15 0.15 0.14 0.13 


2.3 


3.1 


1.4 


2.8 


1.04 


yyny 


+9.87c-04 


± 


1 


7e-04 


1 


J0328+2552 


03:28:44.34 


25:52:08.4 


3.5 


0.16 0.17 0.19 0.17 


1.8 


3.5 


5.0 


1.5 


0.89 


yyyy 


+6.18c-03 


± 


1 


Oc+03 


3 


J0331+2403 


03:31:47.47 


24:03:10.2 


2.4 


0.12 0.12 0.12 0.12 


1.1 


1.2 


2.4 


1.4 


1.22 


nnnn 


+4.71e-04 


± 


7 


2e-05 


2 



Table 4 — Continued 



Name RA(J2000) Dcc(J2000) Ha R SV(Jy) Mod index % %err Var D(2d) tim 



J0338+3818 


03 


:38: 


:27 


.99 


38 


:18: 


56 


,5 


6.2 





07 





.07 


0. 


08 





.10 


1.6 


2.6 


3.3 


1.2 


1 


69 


nnyn 


+1.32e-03 ± 


-1.0e+00 


3 


J0342+3859 


03: 


:42: 


:16 


26 


38: 


:59: 


:06 


,2 


8.1 


0, 


10 





.11 





12 





.13 


3.1 


3.2 


4.5 


3.1 


1 


.24 


yyyy 


+3.83e-03 ± 1.5e-04 


3 


J0343+3622 


03 


:43: 


:28 


95 


36 


:22: 


:12 


,4 


7.8 





75 





.51 





56 





.73 


8.3 


4.2 


2.2 


4.3 


0. 


57 


yyyy 


+7.91c-03 ± 4.0e-04 


3 


J0344+6827 


03: 


:44 


:41 


44 


68' 


:27 


:47. 


,8 


2.5 


0, 


09 





.10 





12 





11 


5.5 


1.5 


3.0 


4.2 


1 


,31 


ynyy 


+5.01c-03 ± 2.5e-04 


3 


J0345+1453 


03: 


:45: 


:06 


41 


14: 


:53 


:49 


,5 


2.2 





33 





.32 





33 





.35 


1.3 


2.0 


1.8 


0.7 





63 


ynyn 


+5.04e-04 ± 3.3e-05 


1 


J0349+4609 


03: 


:49: 


:18 


74 


46: 


:09 


59 


,6 


2.8 


1. 


00 





.97 





91 


0. 


.00 


0.2 


1.0 


0.6 


0.0 


0. 


51 


nnnn 


+1.43e-04 ± 2.4c-05 





J0351+1308 


03' 


:51 


:26 


23 


13: 


:08: 


:05 





3.1 


0, 


10 





.10 


0. 


08 





.10 


2.1 


2.6 


7.8 


6.0 


1 


,48 


nyyy 


+l.llc-02 ± 4.1c-04 


3 


J0354+4643 


03: 


:54: 


:30 


.01 


46: 


.43 


:18 


,7 


2.5 


0. 


,79 





.80 





79 





.00 


1.4 


2.0 


0.5 


0.0 


0. 


,52 


yynn 


+5.01c-04 ± 2.0c-05 


2 


J0355+2005 


03 


:55: 


:19 


85 


20 


05 


:08 





1.9 





,07 





.07 


0. 


06 





.06 


1.3 


2.4 


1.9 


1.6 


2. 


,11 


nnnn 


-1.22c-04 ± 4.1e-05 





J0357+0924 


03: 


:57. 


:46 


63 


09 


:24: 


10 


3 


18.0 





07 





.05 


0. 


05 





.05 


2.9 


5.4 


1.6 


4.0 


2. 


35 


nynn 


+2.49c-03 


± 


3.8e-04 


2 


J0358+3850 


03: 


:58 


:44 


43 


38' 


:50 


:21 


,6 


12.3 





12 





.12 





11 





.11 


1.8 


4.0 


2.3 


1.8 


1 


,24 


nyyn 


+2.97e-03 


± 


1.0e+03 


3 


J0359+6005 


03: 


:59: 


:02 


.64 


60 


:05: 


:22 





4.4 


0. 


,85 





.84 


0. 


78 





.74 


0.6 


1.4 


0.9 


0.8 


0. 


53 


nyny 


+3.18e-04 


± 


-1.0e+00 





J0403+2600 


04: 


:03 


05 


.58 


26: 


:00 


01 


,5 


4.7 


0. 


,86 





.84 





83 


0. 


.83 


0.1 


1.3 


1.4 


0.7 


0. 


52 


nnyn 


+4.20e-04 


± 


2.4c-05 


3 


J0406+2511 


04: 


:06 


:25 


.99 


25: 


:11: 


:31 


,2 


4.5 


0. 


11 





.09 





10 





.11 


2.3 


4.5 


4.5 


8.8 


1 


,39 


nyyy 


+8.63e-03 


± 


4.9c-04 


2 


J0409+1217 


04: 


:09: 


:22 


00 


12 


:17. 


:39 


,8 


7.6 





,73 





.65 





59 





.45 


1.2 


4.1 


4.4 


1.9 


0. 


.55 


yyyy 


+3.52e-03 


± 


2.0c-04 


2 


J0411+0843 


04: 


:11: 


:33 


.85 


08 


:43: 


:11. 


,4 


10.7 





13 





.10 





10 


0. 


.09 


4.7 


5.1 


5.4 


4.0 


1 


.32 


yyyy 


+5.19c-03 


± 


4.0c-04 


1 


J0414+3418 


04: 


:14 


:37 


.25 


34: 


18 


:51 


,1 


13.3 


1 


36 


1 


.37 


1. 


32 


1 


.36 


0.7 


0.2 


0.9 


0.2 


0. 


.51 


nnnn 


+1.49e-04 


± 


1.8c-05 





J0419+3955 


04: 


:19: 


:22 


.55 


39: 


:55: 


:28 


,9 


3.9 


0. 


34 





.37 


0. 


39 





.41 


3.2 


4.2 


2.2 


2.4 


0. 


.60 


yyyy 


+2.54e-03 


± 


l.lc-04 


2 


J0422+0219 


04: 


:22: 


:52 


.21 


02: 


:19 


:26 


9 


4.6 


1 


,09 


1 


.21 


1 


27 


1 


.27 


1.1 


0.2 


1.4 


0.6 


0. 


51 


ynyn 


+3.08e-04 


± 


2.7c-05 





J0423+4150 


04: 


:23: 


:56 


01 


41: 


:50: 


:02 


,7 


2.8 


1 


,70 


1 


.68 


1 


71 





.00 


0.1 


1.0 


1.2 


0.0 


0. 


.51 


nyyn 


+3.09e-04 


± 


-l.Oc+00 





J0427+4133 


04: 


:27: 


:46 


,04 


41: 


33 


:01 





2.8 





,79 





.83 





78 





.00 


1.6 


1.1 


1.0 


0.0 


0. 


.53 


yyny 


+4.17c-04 


± 


2.9e-05 


2 


J0432+3131 


04: 


:32: 


06 


44 


31: 


31 


:13. 


,3 


3.7 


0. 


,12 





.11 





11 


0. 


.12 


0.9 


2.6 


2.4 


2.0 


1 


24 


nyyn 


+1.22e-03 


± 


1.3e-04 


3 


J0443+3441 


04: 


:43 


:31 


.63 


34: 


:41: 


:06 


,6 


2.9 





,87 





.84 





.83 





.81 


0.6 


0.8 


0.9 


0.7 





.52 


nnnn 


+1.18e-04 


± 


9.8c-06 





J0449+1121 


04: 


:49: 


:07 


.67 


11: 


:21 


:28 


,6 


9.8 


2 


.31 


2 


.25 


1 


90 


1 


.89 


1.4 


1.0 


1.3 


0.1 





.51 


yyyn 


+5.37c-04 


± 


l.Oc+03 


3 


J0453+0128 


04: 


53 


02 


.23 


01 


:28: 


:35 


6 


8.4 


0. 


.17 





.18 





17 


0. 


.18 


5.9 6.1 13.3 


, 2.7 


0. 


90 


yyyy 


+2.45e-02 


± 


1.5e-03 


2 


J0459+0229 


04: 


59 


:52 


.05 


02 


:29 


31 


,1 


7.4 


1 


.21 


1 


.21 


1 


22 


1 


.13 


0.4 0.7 1.2 


0.8 


0. 


51 


nnnn 


+1.41e-04 


± 


1.2e-05 





J0502+1338 


05: 


:02 


:33 


21 


13 


38: 


:10 


,9 


8.6 





,54 





.69 





53 





.53 


5.0 4.8 


4.2 


5.7 


0. 


.55 


yyyy 


+9.64e-03 


± 


2.6e-04 


3 


J0503+0203 


05: 


03 


21 


19 


02 


:03: 


:04 


,6 


9.9 


2. 


00 


1 


.95 


1 


94 


1 


.88 


0.1 0.8 


1.0 


0.8 


0. 


.50 


nnnn 


+1.49e-04 


± 


1.8e-05 





J0505+0459 


05 


05: 


:23 


,18 


04: 


59 


:42 


,7 


12.3 





,79 





.76 


0. 


86 


0. 


.79 


1.3 1.3 1.1 


1.0 


0. 


52 


nnnn 


+3.12e-04 


± 


1.4c-05 





J0509+0541 


05: 


:09: 


:25 


,96 


05: 


:41: 


:35 


3 


24.5 


0. 


52 





,49 





56 





.69 


3.1 3.8 6.1 


6.9 





.55 


yyyy 


+1.38e-02 


± 


4.8c-04 


3 


J0510+1800 


05: 


:10: 


:02 


36 


18' 


00 


:41 


,5 


22.1 


0. 


68 





,71 





.58 





.63 


1.8 1.8 1.6 


1.0 


0. 


.55 


yyyn 


+6.60e-04 


± 


6.8c-05 


1 


J0512+2037 


05 


12 


:39 


07 


20 


:37: 


:42 


,7 


6.1 


0. 


,12 





.12 


0. 


11 


0. 


.12 


1.4 4.0 4.1 


1.6 


1 


24 


nyyn 


+3.14e-03 


± 


2.0c-04 


3 


J0530+1331 


05: 


.30 


:56. 


41 


13: 


:31: 


:55 


,1 


70.5 


2. 


,73 


2 


.42 


2. 


86 


3 


.44 


1.1 6.9 0.6 


1.8 


0. 


.51 


nyny 


+4.58e-03 


± 


2.2e-04 


3 


J0534+1047 


05: 


:34: 


:40 


,73 


10' 


.47 


:16 


,5 


147.0 


0. 


14 





.14 





15 





.16 


1.2 0.9 0.9 


0.8 


1 


.01 


nnnn 


+1.52e-05 


± 


1.7c-05 





J0536+5038 


05: 


36: 


:20 


.23 


50 


:38: 


26 


,2 


1.7 


0. 


.94 





.93 





92 


0. 


.91 


1.0 0.7 0.8 


0.3 


0. 


52 


nnnn 


+1.07e-04 


± 


8.0e-06 






Table 4 — Continued 



Name RA(J2000) Dcc(J2000) Ha R S v (Jy) Mod index % %err Var D(2d) tim 



J0539+1433 


05: 


39: 


12 


.36 


14: 


33 


:45. 


5 


30.3 


1. 


04 


1. 


19 


1 


.20 


1 


23 





.1 


1 


.6 


1 


.8 





.2 





.51 


nyyn 


+7.86e-04 


± 


8.9e-05 





J0541+5312 


05: 


41: 


16 


.17 


53: 


12 


:24. 


.8 


1.4 


0. 


43 


0. 


42 





.40 





.41 





.8 





.7 


1 


.1 





.9 





.59 


nnnn 


+1.23c-04 


± l.lc-05 





J0544+4541 


05: 


44: 


01 


.16 


45: 


41 


:02. 


7 


2.1 


0. 


15 





16 





.15 





15 


2 


.5 


1 


.3 


1 


.4 





.9 





.98 


ynnn 


+1.26G-04 


± 1.2c-04 





J0552+0313 


05: 


52: 


50 


.10 


03: 


13 


:27. 


2 


18.4 


0. 


62 


0. 


57 





.56 





.66 





9 


1 


.0 





.9 





.3 





.55 


nnnn 


+1.10e-04 


± 1.5c-05 





J0553+1439 


05: 


53: 


38 


.03 


14: 


39 


:29. 


.7 


9.7 


0. 


14 


0. 


14 





.14 





.13 





6 


2 


.1 


1 


.5 


1 


.1 


1 


.09 


nnnn 


+2.62e-04 


± 4.1c-05 





J0558+0044 


05: 


58: 


.08 


.14 


00: 


14 


:48. 


1 


24.2 


0. 


16 





16 





.17 





17 





.9 


1 


.2 


1 


.1 





.7 





.93 


nnnn 


+8.69e-05 


± 2.5e-05 





J0605+4030 


06: 


05: 


50 


.85 


40: 


30 


:08. 





1.8 


1. 


03 


1. 


06 





.95 





.92 


1 


.3 


1 


.8 


1 


.5 


1 


.1 





.52 


yyyy 


+7.02c-04 


± 2.2e-05 


3 


J0605+5753 


06: 


05: 


42 


.22 


57: 


53 


:16. 


3 


1.3 


0. 


24 


0. 


23 





.22 


0. 


21 





.7 





.9 


1 


.1 


1 


.0 





.76 


nnnn 


+1.49e-04 


± 2.2e-05 





J0606+4401 


06: 


06: 


50 


,21 


44: 


01: 


:40. 


8 


1.6 


0. 


08 





09 





.10 





.09 


1 


.6 


1 


.9 


1 


.8 


1 


.2 


1 


.49 


nnnn 


+3.82e-04 


± 7.4e-05 





J0608+4508 


06: 


08: 


03 


.43 


45: 


08 


;25. 





1.9 


0. 


10 





10 





.10 





10 





.8 


1 


.4 


1 


.4 


1 


.7 


1 


.41 


nnnn 


-9.46e-05 ± 3.5e-05 





J0610+7801 


06: 


10: 


24 


.27 


78: 


01 


:36. 


1 


1.9 


0. 


08 


0. 


08 





.09 





.10 


2 


1 


3 


.5 


1 


.8 


2 


.8 


1 


,61 


nyyn 


+4.19e-03 


± 


1.0e+03 


3 


J0614+6046 


06: 


14: 


23 


.86 


60: 


46 


;21. 


.7 


1.7 


0. 


71 





71 





.72 


0. 


71 


0. 


.8 





.3 





.8 





.8 





.53 


nnnn 


+9.33c-05 


± 


1.0e-05 





J0618+4207 


06: 


18: 


.08 


.61 


42: 


07 


:59. 


8 


1.3 


0. 


16 





15 





.12 





12 


1 


.3 


1 


.2 


2 


.6 


1 


.9 


1 


.07 


ynyn 


+7.33c-04 


± 


5.0e-05 


2 


J0624+3856 


06: 


21: 


19 


.02 


38: 


56 


:48. 


.7 


1.7 


0. 


55 


0. 


55 





.54 


0. 


.60 


1 


.1 


1 


.1 


1 


.0 


1 


.4 





.55 


nnnn 


+3.74e-04 


± 


1.9e-05 





J0639+7324 


06: 


39: 


21 


.96 


73: 


24 


:58. 





1.8 


0. 


78 





.79 





.86 





.84 


1 


.5 


3 


.4 


1 


.3 





.8 





.52 


yynn 


+1.39e-03 


± 


4.5e-05 


3 


J0643+1238 


06: 


43: 


59 


.85 


12: 


38 


;18. 





8.6 


0. 


09 





08 





.10 





.09 


3 


.2 


1 


.4 


3 


.6 


6 


.1 


1 


.47 


yyyy 


+7.71o-03 


± 


1.0c+03 


3 


J0644+2911 


06: 


14: 


11 


.82 


29: 


11 


:04. 





3.5 


0. 


15 


0. 


15 





.15 





.16 


3 


.0 


2 


.0 


2 


.9 


2 


.1 





.99 


ynyn 


+2.09e-03 


± 


1.2e-04 


3 


J0644+3914 


06: 


14: 


53 


.71 


39: 


11 


:47. 


.5 


1.1 





59 





58 





.63 





.67 





6 


1 


.2 


1 


.3 


1 


.1 





.54 


nnyn 


+3.46c-04 


± 


1.6e-05 





J0646+4451 


06: 


46: 


32 


.02 


44: 


51 


:16. 


5 


1.0 


2. 


67 


2. 


69 


2 


.82 


2. 


.90 


0. 


.3 





.5 





.5 





.1 





.50 


nnnn 


+3.76e-06 


± 


1.0e-05 





J0650+6001 


06: 


50: 


31 


.25 


60: 


01 


:44 


5 


1.3 


1. 


16 


1. 


19 


1 


.11 


1 


.10 





.2 


1 


.3 


1 


.0 


1 


.3 





.51 


nyny 


+4.33c-04 


± 


2.1e-05 


3 


J0653+3705 


06: 


53: 


58 


.28 


37: 


05 


:40. 


6 


2.0 


0. 


56 


0. 


57 





.57 


0. 


.60 





.8 


1 


.3 





.7 





.8 





.55 


nnnn 


+1.90e-04 


± 


1.4e-05 





J0657+3858 


06: 


57: 


.51 


.46 


38: 


58 


;02. 


5 


1.1 


0. 


07 





.08 





.07 





08 


2 


.5 


3 


.4 


3 


.5 


4 


.7 


1 


.80 


nnny 


+3.00e-03 


± 


1.7e-04 


2 


J0659+0813 


06: 


59: 


17 


.99 


08: 


13 


:30. 


9 


7.6 


0. 


84 





83 





.85 





.77 





.6 


1 


.1 





.9 





.5 





.52 


nnnn 


+1.55c-04 


± 


1.5e-05 





J0700+1709 


07: 


00: 


01 


.52 


17: 


09 


:21. 


.7 


3.8 


1. 


15 


1. 


09 





.98 





.90 


1 


.0 


1 


.4 


2 


.4 





.5 





.53 


nyyn 


+6.78c-04 


± 


7.2e-05 


2 


J0704+4502 


07: 


04: 


50 


.96 


45. 


02 


:41. 


6 


1.0 





09 





09 





.09 





09 


2. 


.1 


2 


.2 


1 


.8 


1 


.9 


1 


.55 


nnnn 


+3.92c-04 


± 


6.0e-05 





J0708+3455 


07: 


08: 


24 


.44 


34: 


55 


:42. 


1 


2.0 


0. 


15 





15 





.15 


0. 


.17 


1. 


.1 


2 


.5 


2 


.7 


1 


.5 





.97 


nyyn 


+1.33c-03 


± 


1.3e-04 


3 


J0709+3422 


07: 


09: 


:55 


.43 


34: 


22 


:56. 


.7 


2.4 


0. 


10 


0. 


10 





.10 





10 


2 


.2 


2 


.3 


1 


.4 


1 


.6 


1 


.35 


nnnn 


+2.53c-04 


± 


3.0e-05 





J0713+4349 


07: 


13: 


.38 


.16 


43: 


49: 


:17. 


1 


1.2 


1. 


56 


1 


55 


1 


.54 


1 


.54 





.1 





.8 





.1 





.2 





.51 


nnnn 


+1.29e-04 


± 


2.5e-05 





J0713+5053 


07: 


13: 


12 


.89 


50; 


53 


■A3. 


9 


0.5 


0. 


15 





16 





.16 





16 


1 


.5 


1 


.5 


1 


.6 





.9 





.96 


nnnn 


+2.26e-04 


± 


9.8e-05 





J0714+6814 


07: 


14: 


58 


.37 


68: 


11 


:00. 


.5 


1.2 


0. 


.10 





10 





.10 





10 


1 


.8 


1 


.5 


3 


.3 


3 


.5 


1 


.42 


nnyy 


+1.43e-03 


± 


l.le-04 


1 


J0717+5231 


07: 


17: 


29 


.94 


52; 


31 


:04. 


.8 


0.5 


0. 


.14 


0. 


14 





.14 





13 


1 


.0 


1 


.0 


1 


.0 


1 


.0 


1 


.08 


nnnn 


-3.06e-05 ± 1.2e-05 





J0725+1425 


07: 


25: 


16 


.80 


14: 


25 


:13. 


7 


4.1 


0. 


95 


1. 


09 


1 


.12 





.98 


2 


.6 


2 


.9 


1 


.6 


1 


.0 





.51 


yyyn 


+3.97c-03 ± 1.0c+03 


3 


J0726+4124 


07: 


26: 


■22 


.42 


41: 


24 


:43. 


6 


1.1 


0. 


16 


0. 


16 





.14 





13 


1 


.1 


2 


.2 


1 


.6 


1 


.5 


1 


.00 


nnnn 


+5.12c-04 ± 4.7e-05 


1 


J0726+6125 


07: 


26: 


:51 


.67 


61: 


25 


:13. 


6 


1.0 


0. 


16 


0. 


14 





.12 


0. 


.09 


2 


.5 


2 


.1 


1 


.1 


3 


.9 


1 


.17 


ynyy 


+2.82c-03 ± 1.2e-04 


2 


J0726+7911 


07: 


26: 


11 


.74 


79: 


11 


:31. 





1.0 


0. 


.85 


0. 


.88 





.85 


0. 


86 


1 


.5 





.7 


1 


.4 





.8 





.52 


nnyn 


+3.28c-04 ± 1.5e-05 






Table 4 — Continued 



Name RA(J2000) Dcc(J2000) Ha R S v (Jy) Mod index % %err Var D(2d) tim 



J0728+5701 


07: 


28 


:49. 


.63 


57 


:01 


:24 


.3 





9 


0. 


14 


0. 


46 





.49 





.50 


2.0 


0.8 


1.1 


2.2 





.57 


ynny 


+6.69e-04 


± 


3.6e-05 


2 


J0734+6528 


07: 


34 


:18 


.03 


65 


:28 


:02 


.6 


1 





0. 


10 


0. 


.09 





.10 


0. 


.10 


1.7 


1.8 


1.5 


3.0 


1 


.44 


nnny 


+7.16c-04 


± 


1.2o-04 


2 


J0738+1742 


07: 


.38 


:07 


.39 


17 


:42 


:18 


.9 


4. 





1. 


62 


1 


60 


1 


.70 


1 


.71 


0.1 


1.3 


1.1 


0.1 





.50 


nnyn 


+2.81c-04 


± 


4.3e-05 





J0739+0137 


07: 


39 


:18 


.03 


01 


:37 


:04 


.6 


1 


.8 


1. 


43 


1 


44 


1 


56 


1 


61 


1.4 


0.3 


5.5 


1.9 





.50 


ynyy 


+6.84e-03 


± 


9.0e-04 


3 


J0739+1739 


07: 


39 


:28 


.55 


17 


:39 


:27 


.1 


4. 


.0 


0. 


08 





08 





,08 


0. 


.08 


1.3 


2.1 


2.1 


4.0 


1 


.77 


nnny 


+1.87c-03 


± 


1.0c+03 


3 


J0739+7527 


07: 


39 


:13 


.19 


75 


:27 


:47 


.7 


1 


.3 


0. 


21 





21 





.21 





.21 


1.0 


1.0 


1.0 


0.9 





.79 


nnnn 


+9.16e-05 


± 


8.7c-06 





J0741+2557 


07: 


:41 


:29 


.74 


25 


:57 


:32 


.2 


2 


.8 


0. 


07 


0. 


07 





.07 





.07 


1.7 


2.3 


2.3 


1.7 


1 


.96 


nnnn 


+4.82c-04 


± 


-l.Oc+00 


3 


J0741+3112 


07: 


41 


:10 


.70 


31 


:12 


:00 


.2 


2 


3 


2. 


70 


2. 


49 


2 


.25 


2. 


.07 


0.4 


0.2 


1.5 


0.4 





.52 


nnyn 


+4.75c-04 


± 


3.4e-05 





J0745+1011 


07: 


45 


:33 


.05 


10 


:11 


:12 


.6 


1 


.4 


3. 


40 


3 


39 


3 


.40 


3 


.36 


0.1 


1.1 


1.1 


0.1 





.50 


nnyn 


+3.37c-04 


± 


6.2e-05 





J0748+2400 


07: 


48 


:36 


11 


24: 


:00 


:24 


.1 


4. 





0. 


57 





62 





.60 





62 


1.2 


1.1 


1.8 


1.4 





.54 


nnyn 


+4.92c-04 


± 


2.6e-05 


1 


J0749+0837 


07: 


49 


:58 


.53 


08 


:37 


:30 


.7 


1 


3 


0. 


11 





11 





11 





11 


1.2 


2.7 


2.5 


1.8 


1 


.28 


nyyn 


+1.29c-03 


± 


7.3e-05 


3 


J0750+1231 


07: 


50 


:52 


.04 


12 


:31 


:04 


.8 


1 


.8 


1. 


37 


1 


38 


1 


.33 


1 


.34 


0.7 


2.9 


2.1 


2.3 





.51 


nyyy 


+1.30e-03 


± 


1.0e-04 


2 


J0753+5352 


07: 


53 


:01 


.38 


53 


:52 


:59 


.6 





.8 


0. 


71 





96 





.92 


1 


.00 


1.0 


1.5 


2.4 


0.9 





.52 


nnyn 


+7.01c-04 


± 


3.9e-05 


3 


J0757+0956 


07: 


■57 


:06 


.61 


09 


:56 


:34 


.8 


1 


.5 


1. 


.48 


1 


34 


1 


36 


1. 


.57 


5.8 


2.4 


3.0 


2.3 





.51 


yyyy 


+6.55e-03 


± 


3.3e-04 


3 


J0757+6110 


07: 


:57 


:44 


.69 


61 


:10 


:32 


.7 


1 





0. 


11 


0. 


11 





.11 





11 


1.6 


1.6 


1.7 


1.4 


1. 


.31 


nnnn 


+2.69c-04 


± 


4.7e-05 





J0758+0827 


07: 


.58 


:28 


.04 


08 


:27 


:09 


.0 


1 


3 


0. 


10 





10 





.11 





.11 


4.5 


3.5 


8.5 


3.4 


1. 


.33 


yyyy 


+1.23e-02 


± 


5.4e-04 


2 


J0800+4854 


08: 


00 


:34 


.37 


48 


:54 


:24 


.0 





,5 


0. 


11 





12 





11 


0. 


11 


4.0 4.7 1.7 11.4 


1 


,26 


yyny 


+8.55e-03 


± 


8.7e-04 


1 


J0801+6639 


08: 


01 


:36 


.38 


66 


:39 


:09 


.9 





.7 


0. 


10 





11 





11 


0. 


11 


1.5 


2.8 


2.5 


2.3 


1 


.30 


nyyy 


+1.35c-03 


± 


l.le-04 


3 


J0802+1809 


08: 


:02 


:48 


.03 


18 


:09 


:49. 


.2 


3 


,9 


0. 


53 


0. 


52 





.54 





.54 


1.6 


0.9 


1.1 


1.7 





.55 


ynny 


+4.26e-04 


± 


3.2e-05 


2 


J0804+1012 


08: 


04 


:07 


.58 


10 


:12 


:13 


.2 


1 


.5 


0. 


08 


0. 


08 





.08 





08 


4.9 


5.0 


2.8 


3.3 


1. 


.69 


yyyy 


+3.78c-03 


± 


4.2e-04 


1 


J0805+6144 


08: 


05 


:18 


.17 


61 


:44 


:23 


.7 


0. 


,8 


1. 


06 


1 


17 


1 


.22 


1. 


.27 


1.0 


1.1 


0.6 


0.2 





.51 


nnnn 


+9.93e-05 


± 


1.3e-05 





J0807+1352 


08: 


07 


:38 


.50 


13 


:52 


:17 


.3 


1 


,7 


0. 


18 


0. 


18 





.17 





17 


1.6 


1.7 


1.7 


0.9 





.90 


ynnn 


+5.00c-04 


± 


4.1e-05 


2 


J0808+0514 


08: 


08 


:38 


.84 


05 


:14 


:39 


.9 


1 


.0 


0. 


11 


0. 


10 





10 


0. 


.13 


1.8 


2.5 


2.8 


1.6 


1 


.28 


nyyn 


+1.02c-03 


± 


5.7e-05 


1 


J0808+4950 


08: 


08 


:39 


.66 


49 


:50 


:36 


.5 


0. 


,6 


0. 


75 





67 





.64 


0. 


.60 


2.6 


1.2 


0.9 


1.0 





,54 


ynnn 


+4.64e-04 


± 


5.1e-05 


1 


J0810+1010 


08: 


10 


:26 


.42 


10 


:10 


:41 


.0 


1 


.3 


0. 


13 





,13 





.13 





12 


1.0 


1.2 


1.2 


0.7 


1 


14 


nnnn 


-1.90e-05 ± 1.9e-05 





J0811+0146 


08: 


11 


:26 


.70 


01 


:46 


:52 


.2 


1 


,7 


0. 


90 





.85 





.82 


0. 


.86 


0.9 


2.0 


3.7 


3.3 





.53 


nyyy 


+2.91e-03 


± 2.2c-04 


3 


J0818+4222 


08: 


18 


:16 


.00 


42 


:22 


:45 


.4 


1 


.2 


0. 


96 





.89 





.79 


0. 


.87 


1.3 


3.5 


1.5 


2.9 





.53 


nyyy 


+1.59c-03 


± 1.8c-04 


1 


J0821+3107 


08: 


21 


:07 


.61 


31 


:07 


:51 


.1 


2 


3 


0. 


12 





.12 





.11 





11 


1.3 


1.9 


1.8 


1.1 


1 


.26 


nnnn 


+4.96c-04 


± 7.9e-05 


3 


J0823+4138 


08: 


23 


:55 


.81 


41: 


:38 


:10 


.1 


1 


,2 


0. 


13 





,14 





.14 





.14 


1.0 


2.5 


1.7 


1.5 


1. 


.07 


nnnn 


+4.14e-04 


± 3.7c-05 


1 


J0824+5552 


08: 


24 


:47 


.23 


55 


:52 


:42 


.6 





,8 


1. 


21 


1 


,16 


1 


.08 


1. 


03 


0.8 


1.0 


0.2 


0.6 





.52 


nnnn 


+9.78c-05 


± l.lc-05 





J0825+0309 


08: 


25 


:50 


.33 


03 


:09 


:24 


.5 





.9 


0. 


96 





90 


1 


30 


1 


26 


1.2 


2.1 


2.3 


1.4 





.51 


nyyy 


+1.31e-03 


± 1.0c-04 


2 


J0825+0831 


08: 


25 


:04 


.77 


08 


:31 


:11 


.0 


2 


.2 


0. 


14 





,14 





.14 





.14 


1.4 


1.5 


1.5 


0.9 


1 


.07 


nnnn 


+ 1.79e-04 


± 3.7c-05 





J0829+4018 


08: 


29 


:35 


.57 


40 


:18 


:59 


.2 


1 


.1 


0. 


08 





.09 





.10 





11 


9.2 11.3 11.5 


. 17.5 


1. 


.47 


yyyy 


+4.75c-02 


± 5.0c-03 


1 


J0830+2410 


08: 


30 


:52 


.08 


24 


:10 


:59 


.8 


2 


.9 


0. 


98 


1 


.03 


1 


.28 


1. 


34 


2.1 


1.2 1.4 1.0 





.51 


ynyn 


+4.75c-04 


± 4.1e-05 


1 


J0831+0429 


08: 


31 


:48 


.87 


04 


:29 


:39 


.0 


1 


4 


1. 


06 





.96 


1 


.13 


1 


.16 


0.2 1.3 1.0 0.7 





.51 


nnyn 


+1.62e-04 


± 1.4e-05 






Table 4 — Continued 



Name RA(J2000) Dcc(J2000) Ha R &,(Jy) Mod index % %err Var £>(2d) tim 



J0839+0319 


08 


:39 


:49 


.19 


03 


:19 


:53 


.8 


1 


.9 





.63 





.61 





.60 





.58 


1 


.3 





.8 


1 


.1 


0.7 





55 


ynyn 


+2 


62e-04 


± 


1.7e-05 





J0839+2002 


08 


:39 


:10 


.89 


20 


:02 


:07 


.3 


2. 


6 





.11 





11 





.10 





.10 


1 


.1 


2 


.9 


1 


.6 


1.1 


1 


33 


nynn 


+7. 


88c-04 


± 


l.Oc+03 


3 


J0839+4301 


08 


:39: 


:52 


.37 


43 


:01 


:48 


.8 





.8 





.10 





10 





.10 





.10 


1 


.6 


2 





1 


.4 


1.5 


1 


.41 


nnnn 


+1. 


13e-04 


± 


2.6e-05 





J0842+1835 


08 


:42 


:05 


.09 


18 


:35 


:10 


.9 


2 


.0 


1 


.08 


1. 


07 


1 


.04 


1 


.02 





.1 


1 


.2 


1 


.2 


0.1 





.52 


nnyn 


+4. 


89e-04 


± 


1.9e-05 





J0850+5159 


08 


:50 


:42 


.24 


51 


:59 


:11 


.6 





.7 





11 


0. 


11 





.11 





.10 


3. 


,5 


1 


.4 


6 


.3 


5.3 


1 


.31 


ynyy 


+8. 


62e-03 


± 


-l.Oc+00 


3 


J0851+4847 


08 


:51 


:36 


.51 


48 


:47: 


:30 


.7 





.8 





.06 





06 





.07 





.08 


2 


,7 


2 


,5 


2 


.9 


1.8 


1 


92 


nnyn 


+1. 


05e-03 


± 


l.le-04 


2 


J0854+2006 


08 


:54 


:48 


.87 


20 


:06 


:30 


6 


1 


.3 


2. 


.28 


2. 


67 


2 


.17 


1 


.93 


1. 


.4 


0. 


.1 





.2 


0.9 





.51 


ynnn 


+3. 


.Olc-04 


± 


4.0c-05 





J0854+8034 


08 


:54 


:48 


.58 


80 


:34 


:22 


.3 


1 


6 





.15 





16 





.17 





.18 


2 


,1 


1 


.5 


1 


.7 


2.5 





92 


yyyy 


+7. 


81c-04 


± 


4.3e-05 


1 


J0856+7146 


08 


:56 


;5d 


.86 


71 


:46 


:23 


.8 





.9 





.20 





17 





.14 





.14 


3 


4 


2 


3 


4 


.0 


5.1 





.94 


yyyy 


+4. 


56e-03 


± 


2.4e-04 


3 


J0902+0443 


09 


:02 


:27 


.16 


04 


:43 


:09 


.5 


2 


.7 





.13 





13 





.14 





.13 


1 


,1 


1. 


.7 


1 


.4 


1.4 


1 


.08 


nnnn 


+4 


71c-04 


± 


5.5e-05 


3 


J0903+6757 


09 


:03 


:53 


.15 


67 


:57: 


:22 


6 





.8 





.60 


0. 


61 





.62 





.62 


1. 


,1 


0. 


.9 


1 


.3 


0.8 





54 


nnyn 


+2. 


77c-04 


± 


2.6e-05 





J0909+0121 


09 


:09 


:10 


.09 


01 


:21: 


:35 


.6 


1 


.4 


1 


.99 


1 


.80 


1 


.85 


1 


.97 


1 





1 


,1 





.1 


0.4 





51 


nnnn 


+1. 


42e-04 


± 


1.7c-05 





J0914+0245 


09 


:14 


:37 


.91 


02 


:45: 


:59 


.2 


2 


.1 





.41 





.41 





.45 





.52 


5 


.2 


1 


,8 


1 


.0 


1.3 





58 


yyyy 


+2. 


03e-03 


± 


1.0c+03 


3 


J0916+0242 


09 


:16 


:41 


.77 


02 


:42: 


:52 


.9 


1 


.9 





.12 





12 





.12 





.14 


5 


.7 


4. 


.4 


2 


.8 


2.2 


1 


17 


yyyy 


+5. 


75e-03 


± 


2.8e-04 


3 


J0920+4441 


09 


:20 


:58 


.45 


44 


:41 


:53 


.9 





.4 


1 


.06 


1 


02 


1 


.03 


1 


.03 





.4 


0. 


9 


1 


.0 


0.5 





.52 


nnnn 


+1. 


12e-04 


± 


1.5e-05 





J0921+6215 


09 


:21: 


:36 


.23 


62 


:15: 


:52 


.1 





.8 


1 


.40 


1 


27 


1 


.17 


1 


.12 


1 





1 


.0 





.8 


0.3 





.51 


nnnn 


+1. 


45e-04 


± 


1.5e-05 





J0922+4209 


09 


:22 


:39 


.58 


42 


:09 


:53 


.3 





.4 





10 


0. 


10 





.10 





.11 


1 


,5 


3 


,4 


1 


.9 


2.1 


1 


.36 


nynn 


+7. 


07c-04 


± 


1.0e-04 


1 


J0923+2018 


09 


:23 


:09 


.07 


20 


:18: 


:49 


A 


1 


.1 





06 


0. 


06 





.07 





.08 


2 


,2 


2 


.3 


3 


.1 


4.5 


1 


.94 


nnny 


+3. 


57c-03 


± 


1.9e-04 


3 


J0929+5013 


09 


:29 


;15 


.44 


50 


:13: 


:35 


.9 





.6 





.51 





.54 





.58 





.45 


6 


.2 


3 


,1 


3 


.4 


5.0 





55 


yyyy 


+4. 


19c-03 


± 


1.7e-04 


1 


J0931+2750 


09 


:31 


:51 


.78 


27 


:50: 


:50 


6 





.8 





.16 





17 





.16 





.16 


1. 


.4 


1 


.8 


3 


.0 


2.1 





95 


nyyy 


+1 


73c-03 


± 


1.2e-04 


3 


J0943+2028 


09 


:43 


:48 


.09 


20 


:28: 


:09 


.9 





6 





.14 





14 





.13 





.14 


1. 


9 


2 


.1 


2 


.3 


1.3 


1 


08 


nynn 


+1. 


Olc-03 


± 


5.9c-05 


3 


J0946+5020 


09 


:46 


:16 


.04 


50 


:20 


:09 


.3 





.5 





.13 





14 





.16 





.16 


8. 


.1 


7. 


.1 


1 


.9 


3.6 


1 


.01 


yyyy 


+8. 


85e-03 


± 


9.0e-04 


1 


J0948+4039 


09 


:48 


:55 


.34 


40 


:39: 


:44 


.5 





.5 


1 


.84 


1 


81 


1 


.84 


1 


.89 


1. 


.0 


1 


3 





.9 


0.7 





50 


yynn 


+1 


80e-04 


± 


2.1c-05 





J0949+5819 


09 


:49 


:39 


.81 


58 


:19: 


:12 


.9 





.7 





.08 


0. 


05 





.03 





.02 


15. 


.1 


5 


,7 


6 


.1 


13.5 


3 


37 


yyny 


+1. 


90e-02 


± 


1.9e-03 


1 


J0952+4941 


09 


:52 


:00 


.66 


49 


:41: 


:25 


.9 





.5 





.15 





15 





.15 





.14 


1 


6 


1 


,1 





.7 


1.0 


1 


02 


nnnn 


+3 


39e-05 


± 


3.2e-05 





J0953+1720 


09 


:53 


:59 


.23 


17 


:20: 


:56 


.6 





.9 





.15 





13 





.13 





.17 


1 


.8 


1 


3 


3 


.9 


2.7 


1 


.03 


ynyy 


+2. 


04e-03 


± 


1.0e-04 


2 


J0956+2515 


09 


:56 


:49 


.87 


25 


:15: 


:16 


.0 





.6 





.97 





93 





.93 





.92 





.1 


1 


.8 


2 


.1 


0.2 





.52 


nyyn 


+1. 


21c-03 


± 


7.7c-05 





J0958+4725 


09 


:58 


:19 


.67 


47. 


:25. 


:07 


.8 





.7 


1 


.67 


1 


.75 


1 


.95 


1 


.77 


2. 


.6 


1 


3 


1 


.3 


1.0 





50 


yyyy 


+7. 


86e-04 


± 


3.4e-05 


2 


J0958+6533 


09 


:58 


:47. 


.24 


65 


:33: 


:54 


.8 


1 


.1 





.70 





62 





.53 





.51 


1 


,1 


1 


9 


1 


.7 


1.4 





,55 


nyyy 


+4. 


80e-04 


± 


5.4e-05 


1 


J1007+1356 


10 


:07 


:41. 


.49 


13 


:56: 


:29 


6 


1 


.2 





.83 





.88 





.85 





.85 





.9 


1 


.3 


1 


.5 


0.7 





52 


nyyn 


+3. 


20e-04 


± 


2.7c-05 





J1008+0621 


10 


:08 


:00 


.81 


06 


:21: 


:21 


.2 





.8 





.72 





83 





.87 





.59 


2. 


.6 


2 


.2 


2 


.5 


2.4 





.53 


yyyy 


+1 


57c-03 


± 


6.5e-05 


1 


J1009+0622 


10 


:09 


:49 


.80 


06 


:22 


:00 


.9 





.8 





.23 


0. 


23 





.23 





.23 





,7 


0. 


,8 


1 


.7 


1.0 





75 


nnyn 


+2. 


92e-04 


± 


2.7c-05 





J1009+1322 


10 


:09 


:00 


.09 


13 


:22 


:54 


.4 


1 


.2 





.07 





07 





.07 





.07 


2. 


,1 


1 


.9 


2 


.8 


2.0 


1 


91 


nnnn 


+7. 


34e-04 


± 


8.3e-05 


3 


J1011+6529 


10 


:11: 


:38 


.18 


65 


:29: 


:21 


.3 





.8 





.12 





12 





.11 





.10 


1. 


.7 


1 


.9 


1 


.9 


2.2 


1 


25 


nnyn 


+5. 


99e-04 


± 


6.1c-05 


1 


J1013+2449 


10 


:13: 


:53 


.42 


24 


:49 


:16 


.4 





.3 





.92 





92 





.93 





.91 





.4 


1. 


6 


2 


.0 


0.8 





52 


nyyn 


+5. 


66e-04 


± 


6.0e-05 


2 
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J1016+2037 


10 


:16 


:44 


.32 


20: 


37 


:47 


.3 


0.8 


0. 


79 





.77 


0. 


77 





78 


0. 


.5 


1 


9 


1 


.6 


1 


1 





.53 


nyyn 


+5 


,22e-04 


± 


3 


4e-05 


2 


J1018+3542 


10 


:18 


:10 


.98 


35: 


42 


:39 


.4 


0.6 


0. 


83 





86 


0. 


90 





96 





9 


0. 


.3 


1 


.5 





,8 





.52 


nnyn 


+2 


90c-04 


± 


2. 


5e-05 





J1021+5114 


10 


:21 


:32 


.23 


51: 


14: 


:33 


.7 


0.3 


0. 


10 





10 


0. 


10 





.11 


1 


.9 


2 


.1 


1 


.9 


1 


,4 


1 


.36 


nnnn 


+3 


,99e-04 


± 


6 


3c-05 





J1022+4126 


10 


:22 


:02 


02 


41: 


26 


:05 


.3 


0.6 





11 


0. 


12 


0. 


13 





14 


2 


.2 


4. 


.2 


2 


.7 


2 


,0 


1 


.16 


nyyn 


+2 


32e-03 


± 


1 


2e-04 


2 


J1024+1912 


10 


:24 


:44 


.81 


19: 


12 


:20 


.4 


0.8 





52 


0. 


51 


0. 


51 


0. 


.56 





.7 


2 


.0 


1 


.2 


1 


.1 





.56 


nynn 


+3 


36e-04 


± 


3 


,7c-05 





J1024+2332 


10 


:24 


:53 


.63 


23: 


32 


:33 


.9 


0.8 


0. 


15 





20 


0. 


17 





.16 


5 


.6 


5. 


.7 


3 


.0 


2 


.2 





.91 


yyyy 


+5 


51e-03 


± 


3 


,4c-04 


2 


J1025+1253 


10 


:25 


:56 


.28 


12: 


53 


:49 


.0 


1.3 


0. 


78 


0. 


78 


0. 


76 


0. 


.69 


1 


.0 


1 


.5 


1 


.4 





.6 





.53 


nyyn 


+4 


09e-04 


± 


2. 


5e-05 


3 


J1026+7032 


10 


:26 


:36 


.68 


70: 


32 


:43 


.0 


1.2 


0. 


12 





12 


0. 


13 


0. 


12 


1 


6 


2. 


.4 


2 


.8 


1 


.9 


1 


.14 


nyyy 


+2 


14e-03 


± 


2 


,2e-04 


3 


J1029+2911 


10 


:29 


:11 


,47 


29: 


11: 


:01 


.9 


0.7 


0. 


09 





08 


0. 


07 


0. 


.07 


2 


.3 


1 


4 


2 


.6 


2 


,7 


1 


.73 


nnnn 


+6 


,74c-04 


± 


9 


5e-05 


1 


J1034+6832 


10 


:3d 


:01 


.11 


68: 


32: 


:27 


.1 


0.7 


0. 


10 





11 


0. 


11 





.10 


1 


.4 


1 


4 


1 


.2 


1 


,3 


1 


.33 


nnnn 


+1 


.89e-05 


± 


3 


3e-05 





J1041+0610 


10 


:41 


:17 


.16 


06: 


10 


:16 


.9 


1.5 


1. 


47 


1. 


40 


1. 


36 


1 


.31 





.5 





.2 


1 


.4 





.2 





.51 


nnnn 


+2 


27c-04 


± 


2 


2e-05 





J1041+4138 


10 


:41 


:16 


.62 


41: 


38. 


:59 


.0 


0.6 


0. 


09 





09 


0. 


10 





11 


1 


.4 


2 


.4 


2 


.5 


2 


.4 


1 


.39 


nnyn 


+ 1 


01c-03 


± 


8. 


le-05 


2 


J1042+3615 


10 


:42: 


:32 


,24 


36: 


15 


:20 


.8 


0.5 


0. 


10 





09 


0. 


09 





09 


2 


.3 


3 


,7 


2 


.1 


1 


,4 


1 


.46 


nnnn 


+ 1 


19c-03 


± 


1 


9e-04 


2 


J1043+2408 


10 


:43 


:09 


,03 


24: 


.08 


:35 


.4 


0.6 


0. 


98 


1. 


02 


0. 


96 





91 





1 





,2 


3 


.5 


1 


.2 





.52 


nnyy 


+2 


52c-03 


± 


1 


,0c+03 





J1044+8054 


10 


:44 


:23 


,06 


80: 


54 


:39 


.4 


1.0 


1. 


09 


1. 


13 


1. 


20 


1. 


.18 


0. 


.8 





.1 





.2 





.2 





.51 


nnnn 


+6 


81e-05 


± 


7 


7c-06 





J1048+7143 


10 


:48 


;27 


.62 


71: 


43 


:35 


.9 


0.5 


1. 


47 


1 


.47 


1. 


30 


1. 


33 


1 


.1 





.9 


1 


.8 





.2 





.51 


nnyn 


+6 


,40e-04 


± 


3 


9e-05 


3 


J1049+1429 


10 


:49 


:46 


.32 


14: 


29 


:38 


.5 


0.9 


0. 


16 


0. 


15 


0. 


15 





.20 


4. 


.3 


3 


,7 


3 


.5 


5 


.9 





.94 


yyyy 


+8 


39e-03 


± 


4 


,7e-04 


3 


J1054+3928 


10 


:54 


:32 


.42 


39: 


28 


:12 


.3 


0.3 


0. 


12 





14 


0. 


.14 





14 


2 


A 


4 


.2 


4 


.3 


2 


.2 


1 


.08 


nyyn 


+3 


60e-03 


± 


3 


6e-04 


1 


J1056+7011 


10 


:56 


:53 


.61 


70: 


11 


:45 


.9 


-0.2 


0. 


53 





.43 


0. 


38 


0. 


.41 


1 


.4 


1 





1 


.1 


2 


.1 





.59 


nnny 


+3 


,97c-04 


± 


3 


4e-05 





J1058+0133 


10 


:58 


:29 


.60 


01: 


33 


:58 


.8 


0.9 


3. 


41 


3 


40 


3. 


14 


3 


10 





1 





6 





.9 





.7 





.50 


nnnn 


+9 


,99e-05 


± 


-1.0e+00 





,11101+0215 


11 


:01 


:37 


.43 


02: 


15 


:12 


.2 


1.0 


0. 


10 





10 


0. 


11 


0. 


12 


2 


.3 


2. 


,8 


2 


.6 


2 


.0 


1 


.30 


nyyy 


+9 


98c-04 


± 


9 


7e-05 


1 


J1104+3812 


11 


:04 


;27 


.31 


38: 


12 


:31 


.7 


0.6 





58 





60 


0. 


65 


0. 


.69 


0. 


.3 


2 


.4 


1 


.5 





.4 





.54 


nyyn 


+7 


26e-04 


± 


5 


,6e-05 





J1105+1754 


11 


:05 


:38 


.81 


17: 


54 


:15 


.6 


0.4 


0. 


14 


0. 


14 


0. 


.14 





.14 


1 


.0 


1 


.9 


1 


.2 


1 


.3 


1 


.05 


nnnn 


+2 


51e-04 


± 


3 


le-05 





J1107+5219 


11 


:07 


:25 


.83 


52: 


19 


:31 


.6 


0.1 


0. 


14 


0. 


14 


0. 


15 





14 


2 


.9 


1 


.8 


2 


.3 


1 


,4 


1 


.04 


ynyn 


+1 


07o-03 


± 


7. 


,7c-05 


2 


J1112+3527 


11 


:12 


:08 


.08 


35: 


27 


:07 


.6 


0.6 





09 





09 


0. 


09 





09 


2 


.1 


2 


.5 


1 


.9 


1 


,6 


1 


.54 


nnnn 


+4 


,40e-04 


± 


6 


le-05 


1 


J1118+2922 


11 


:18 


:57 


.92 


29: 


22 


:13 


.7 


0.5 





09 





09 


0. 


09 





.09 


1 


.3 


2 


.6 


2 


.2 


1 


.7 


1 


.57 


nnnn 


+3. 


53c-04 


± 


4 


,8e-05 





,11120+2108 


11 


:20 


:11 


.65 


21: 


08 


:48 


,1 


0.6 


0. 


09 





09 


0. 


09 





.09 


1 


.1 


3 


.1 


2 


.1 


1 


.7 


1 


.48 


nnnn 


+2 


71e-04 


± 


1 


,7c-04 





J1125+2610 


11 


:25 


:53 


.71 


26: 


10 


:19 


.9 


0.5 


0. 


82 





83 


0. 


83 


0. 


85 





.1 


1 


.8 





.4 





.2 





.52 


nynn 


+7 


08c-04 


± 


9 


,7c-05 





J1130+3815 


11 


:30 


:53 


.28 


38: 


15 


:18 


.5 


0.5 


0. 


.95 





.98 


1. 


03 


1 


.04 


1 


.3 


2 


.0 





.2 


2 


.8 





.52 


nyny 


+2 


18e-03 


± 


1 


,2e-04 


3 


J1132+4455 


11 


:32 


:23 


.76 


44: 


55 


:10 


.3 


0.5 


0. 


13 





13 


0. 


13 





13 


2 


.5 


2 


.8 


1 


.9 


3 


,0 


1 


.13 


nynn 


+1 


.38e-03 


± 


1 


Oe-04 


1 


J1134+2901 


11 


:34 


:14 


.32 


29: 


01 


:21 


.1 


0.4 


0. 


16 


0. 


16 


0. 


16 





16 





.6 


1 


.6 


1 


.3 





.9 





.95 


nnnn 


+1 


26e-04 


± 


2. 


,7e-05 





J1137+3441 


11 


:37 


:09 


.16 


34: 


41 


:55 


.9 


0.5 


0. 


08 





08 


0. 


08 





.08 


3 


.2 


3 


.9 


2 


.1 


1 


,8 


1 


.72 


nynn 


+ 1. 


32e-03 


± 


1 


,3e-04 


1 


J 1141+5953 


11 


:41 


:16 


.06 


59: 


53 


:08 


.7 


0.4 


0. 


.11 





11 


0. 


12 


0. 


12 


1 


3 


2. 


,5 


1 


.7 


1 


.5 


1 


.25 


nnnn 


+4 


,44e-04 


± 


2. 


9e-05 


1 


J1149+3559 


11 


:49 


:33 


.97 


35: 


59 


:08 


,0 


0.5 


0. 


14 





15 


0. 


18 


0. 


.15 


1 


.5 


2 


.9 


2 


.3 


1 


.3 





.96 


nyyn 


+1 


21e-03 


± 


6 


,8e-05 


1 


J1150+2417 


11 


:50: 


:19 


,21 


24: 


17. 


:53 


.8 


0.4 


0. 


.72 





72 


0. 


72 





.78 


0. 


.5 


1 


.0 


1 


.0 





.6 





.53 


nnnn 


+9 


68e-05 


± 


1 


, le-05 
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J1153+8058 


11 


:53 


:12 


.49 


80 


:58: 


29. 


1 


1 


.2 


1. 


47 


1 


56 


1 


.57 


1 


57 


0. 


.2 





.9 


1 


.2 





.6 





.51 


nnyn 


+2.47e-04 ± 6.5e-06 





J1157+0641 


11 


:57 


:00 


.65 


06 


:41: 


12. 


6 





.8 


0. 


18 


0. 


18 





.16 





14 


1 


.7 


1 


.4 


1 


.9 


1 


.5 





.94 


nnyn 


+5.25c-04 ± 3.5e-05 


1 


J1158+2450 


11 


:58 


;25 


.78 


24 


:50: 


17. 


9 





.3 


0. 


80 





.79 





.78 





77 





.3 





.3 





.6 





.4 





.53 


nnnn 


-1.44e-05 ± - 


1.0e+00 





J1159+2914 


11 


:59 


:31 


.83 


29 


:14: 


43. 


.8 





.5 


2. 


89 


2. 


89 


2 


.28 


1. 


.90 


5 


.8 


3 


.1 


4 


.6 


4 


.5 





.52 


yyyy 


+6.27e-03 


± 


5.4e-04 


1 


J1203+1118 


12 


:03 


:12 


.98 


11 


:18: 


17. 


3 





.8 


0. 


17 


0. 


21 





.22 





.23 


1 


1 


1 


.7 


1 


.9 


1 


.1 





.80 


nnnn 


+5.65e-04 


± 


4.8e-05 


2 


J1219+4829 


12 


:19 


;06 


.41 


48 


:29: 


56. 


1 





.5 


0. 


70 





66 





.61 


0. 


59 





.6 


1 


.2 





.9 





.5 





.54 


nnnn 


+1.58e-04 


± 


1.5e-05 





J1221+2813 


12 


:21 


:31 


.69 


28 


:13: 


58. 


.4 





.4 


0. 


74 


0. 


.70 





.66 


0. 


64 


1 


.4 


1 


.0 





.6 





.3 





.54 


ynnn 


+2.13e-04 


± 


2.9e-05 





J1222+0413 


12 


:22 


:22 


.55 


04 


:13: 


15. 


.7 





.2 


0. 


62 


0. 


57 





.59 


0. 


62 





.7 


1 


.2 


1 


.2 





.5 





,55 


nnnn 


+1.66c-04 


± 


1.5e-05 





J1224+0330 


12 


:24 


:52 


.42 


03 


:30: 


50. 


2 





.2 


0. 


82 





.80 





.78 


0. 


76 


1 


.0 





.7 


1 


.3 





.4 





.53 


nnnn 


+1.38c-04 


± 


1.2e-05 





J1229+5522 


12 


:29 


;09 


.28 


55 


:22: 


30. 


5 





.3 


0. 


07 





09 





.11 





12 


2 


.8 


1 


.6 


3 


.5 


2 


.2 


1 


.43 


nnyn 


+2.18e-03 


± 


1.8e-04 


3 


J1242+3751 


12 


:42 


:51 


.37 


37 


:51: 


00. 








.5 


0. 


65 


0. 


66 





.67 





.70 


1. 


.4 





.3 


1 


.2 





.9 





.54 


ynnn 


+3.93e-04 


± 


3.4e-05 





J1247+7046 


12 


:47: 


;07 


.55 


70 


:46: 


45. 


1 





.4 


0. 


13 





14 





.17 





.17 


2 


.6 


2 


.7 


3 


.3 


2 


.0 





.98 


yyyy 


+2.85e-03 


± 


1.0e+03 


3 


J1247+7124 


12 


:47 


:09 


.32 


71 


:24: 


20. 








.4 


0. 


13 





13 





.13 





13 


1 


.4 


1 


.3 


1 


.2 


1 


.5 


1 


.10 


nnnn 


+1.52e-04 


± 


2.5e-05 





J1248+2022 


12 


:48 


:37 


.26 


20 


:22: 


26. 


3 





.7 


0. 


08 


0. 


08 





.08 


0. 


08 


1 


.9 


1 


.3 


3 


.0 


3 


.3 


1 


.74 


nnyn 


+8.79e-04 


± 


1.3e-04 


1 


J1252+1607 


12 


:52 


:36 


.51 


16 


:07: 


06. 


9 





.3 


0. 


11 





11 





.11 





12 


1 


.1 


2 


.2 


1 


.4 


1 


.8 


1 


.26 


nnnn 


+5.70e-04 


± 


6.1e-05 


3 


J1254+0233 


12 


:54 


:45 


.46 


02 


:33: 


.28. 


9 





.4 


0. 


22 





22 





.21 





20 





.8 


1 


.2 


1 


.3 





.8 





.79 


nnnn 


+2.56e-04 


± 


1.0e+03 





J1254+1141 


12 


:54 


:38 


.25 


11 


:41: 


05. 


.8 





.3 


0. 


64 


0. 


65 





.67 





.67 





6 





.8 


1 


.6 


1 


.0 





.54 


nnyn 


+3.40e-04 


± 


1.0e+03 





J1254+1856 


12 


:54 


:33 


.27 


18 


:56: 


01. 


9 





.6 


0. 


13 


0. 


13 





.13 





12 


1 


1 





.8 


2 


.0 


1 


.7 


1 


.14 


nnnn 


+2.57e-04 


± 


4.0e-05 





J1256+1247 


12 


:56 


:34 


.85 


12 


:47: 


34. 


2 





.5 


0. 


.10 





09 





.08 





08 


1. 


.1 


2 


.0 


3 


.2 


1 


.9 


1 


.57 


nyyy 


+1.01o-03 


± 


1.4e-04 


3 


J1300+0828 


13 


:00 


:36 


.44 


08 


:28: 


02. 


.8 





.3 


0. 


12 


0. 


12 





.12 





12 


3 


.8 


2 


.1 


2 


.3 


1 


.2 


1 


.20 


yyyn 


+2.30e-03 


± 


1.0e+03 


3 


J1302+4819 


13 


:02 


:17 


.19 


48 


:19: 


17. 


5 





.5 


0. 


29 


0. 


27 





.24 


0. 


23 


1 


.0 


1 


.8 


1 


.6 





.9 





.72 


nnnn 


+3.54e-04 


± 


3.7e-05 





J1302+5748 


13 


:02 


:52 


.46 


57 


:48; 


37. 


6 





.3 


0. 


95 





89 





.84 


0. 


.91 


1 


.3 





.3 





.9 





.8 





.52 


nnnn 


+1.72e-04 


± 


2.6e-05 





J1309+0716 


13 


:09 


:29 


.18 


07 


:16: 


56. 


3 





.6 


0. 


.10 





10 





.10 





11 


1 


.7 


1 


.6 


1 


.8 


2 


.3 


1 


.37 


nnnn 


+3.58o-04 


± 


4.3e-05 





J1309+1154 


13 


:09 


:33 


.93 


11 


:54: 


:24. 


5 





.4 


1. 


13 


1. 


15 


1 


.10 


1. 


.07 





3 





.2 


1 


.5 





.2 





.51 


nnnn 


+4.12e-04 


± 


5.5e-05 





J1310+3220 


13 


:10 


:28 


.66 


32 


:20: 


43. 


7 





.2 


1. 


63 


1. 


.74 


1 


.81 


1 


74 





.8 





.1 


1 


.2 





.9 





.51 


nnnn 


+2.20e-04 


± 


3.8e-05 





J1310+3233 


13 


:10 


:59 


.40 


32 


:33: 


34. 


4 





.2 


0. 


69 





.72 





.74 


0. 


.79 





.7 


1 


.3 


1 


.3 


1 


.2 





.53 


nynn 


+2.42e-04 


± 


2.9e-05 





J1312+3113 


13 


:12 


:48 


,80 


31 


:13: 


00. 


.5 





.3 


0. 


10 


0. 


10 





.12 





13 


1 


.6 





.9 


1 


.8 


1 


.7 


1 


.28 


nnnn 


+3.26c-04 


± 


6.1e-05 





J1316+0504 


13 


:16 


;22 


.95 


05 


:04: 


52. 


6 





.5 





09 


0. 


09 





.09 





10 


2 


.0 


1 


.5 


2 


.6 


2 


.7 


1 


.48 


nnny 


+1.17c-03 


± 


l.lc-04 


3 


J1316+6927 


13 


:16 


:22 


.98 


69 


:27: 


16. 


6 





.4 


0. 


.11 





11 





.12 





12 


1 


.5 


1 


.3 


1 


.8 


1 


.6 


1 


.24 


nnnn 


+4.19e-04 


± 


5.2e-05 


3 


J1327+2210 


13 


:27 


:00 


.86 


22 


:10: 


50. 


1 





.5 


1. 


18 


0. 


98 





.85 





.82 





1 





.9 


1 


.4 





.4 





.53 


nnnn 


+2.94e-04 


± 


2.5e-05 





J1328+6221 


13 


:28 


:10 


.56 


62 


:21: 


37. 








.5 


0. 


.12 





.12 





.12 





12 


5 


.0 


6 


.0 


8 


.0 


5 


.2 


1 


.17 


yyyy 


+8.36e-03 


± 


5.0e-04 


1 


J1329+3154 


13 


:29 


;52 


.86 


31 


:54: 


11. 








.2 


1. 


47 


1 


32 


1 


.15 


1 


12 





.6 





.9 





.2 





.5 





.51 


nnnn 


+9.31e-05 


± 


1.3e-05 





J1330+4954 


13 


:30 


:29 


.38 


49 


:54; 


27. 


.4 





.5 


0. 


12 





.12 





.12 





11 


2 


.1 


1 


.6 


1 


.5 


1 


.2 


1 


,23 


nnnn 


+2.36c-04 


± 


4.5e-05 





J1350+3034 


13 


:50 


:52 


.73 


30 


:34: 


53. 


5 





.4 


0. 


33 





.32 





.26 





30 


1 


.0 


2 


.0 


1 


.7 


1 


.1 





.67 


nyyn 


+6.99e-04 


± 


4.6e-05 


2 


J1352+3603 


13 


:52 


:00 


.94 


36 


:03: 


51. 


2 





.4 


0. 


.11 


0. 


10 





.10 





10 


1 


.4 


1 


.2 


1 


.7 


1 


.0 


1 


.33 


nnnn 


+1.73e-04 


± 


3.6e-05 






Table 4 — Continued 



Name RA(J2000) Dcc(J2000) Ha R S v (Jy) Mod index % %err Var D(2d) tim 



J1354+6645 


13: 


54 


:23 


.08 


66 


:45 


:25. 


5 


0. 


6 


0. 


12 





13 


0. 


12 





11 


3 


.9 


3 


.9 


3 


.2 


1 


.7 


1 


.19 


yyyn 


+3.26c-03 ± 1.8c-04 


2 


J1355+3024 


13: 


55 


:41 


.11 


30 


:24 


:11. 


.8 


0. 


5 


0. 


11 


0. 


11 


0. 


11 


0. 


11 


1 


.1 


1 


.5 


2 


.0 


1 


.0 


1 


.30 


nnyn 


+1.94e-04 ± 4.6e-05 





J1407+2827 


14: 


07 


:00 


.39 


28 


:27 


:14. 


6 


0. 


4 


2. 


46 


2. 


.47 


2. 


47 


2. 


.48 





.1 





.1 





.7 





.0 





.50 


nnnn 


+5.44c-05 


± 1.8e-05 





J1410+6141 


14: 


10 


:30 


.98 


61 


:41 


:36. 


9 


0. 


4 


0. 


13 


0. 


13 


0. 


14 





14 


1 


3 





.8 


1 


.5 


1 


.0 


1 


.09 


nnnn 


+2.46e-04 ± 3.8c-05 





J1414+1922 


14: 


14 


:29 


.50 


19 


:22 


:18. 


.4 


0. 


7 


0. 


16 





15 


0. 


15 





.15 





6 





.8 


1 


.5 





.8 





.98 


nnnn 


+1.79c-04 ± 1.0c+03 





J1415+1320 


14: 


15 


:58 


.81 


13 


:20 


:23. 


7 


0. 


.4 


0. 


85 





89 


0. 


94 





92 





.7 


1 


.0 


1 


.2 





.8 





.52 


nnyn 


+2.58c-04 ± l.le-05 





J1416+1242 


14: 


16 


:28 


.65 


12 


:42 


:13. 


5 


0. 


.4 


0. 


10 


0. 


10 


0. 


10 





.10 


1 


.3 


1 


.7 





.7 


1 


.3 


1 


.39 


nnnn 


-1.82c-05 ± 4.3e-05 





J1417+3818 


14: 


17 


■AO 


.44 


38 


:18 


:21. 


1 


0. 


3 


0. 


10 


0. 


09 


0. 


12 





12 


2 


.7 


2 


.5 


2 


.6 


2 


.7 


1 


.32 


yyyy 


+1.45e-03 


± 


1.4e-04 


1 


J1419+3821 


14: 


19: 


:46 


.61 


38 


:21: 


:48. 


4 


0. 


3 


0. 


60 





67 


0. 


66 


0. 


73 


1 


.0 





.5 


1 


.1 





.9 





.54 


nnnn 


+2.10o-04 


± 


2.2e-05 





J1419+5423 


14: 


19 


:46 


.59 


54 


:23 


:14. 


7 


0. 


5 


0. 


95 


0. 


89 


0. 


95 





91 


1 


.3 


1 


.2 


1 


.1 





.7 





.52 


nyyn 


+2.16e-04 


± 


1.8e-05 





J1426+3625 


14: 


26 


:37 


.08 


36 


:25 


:09. 


5 


0. 


2 


0. 


41 


0. 


42 


0. 


42 





46 


1 


.0 


1 


.5 


1 


.0 


1 


.3 





.59 


nnnn 


+2.87c-04 


± 


3.7c-05 





J1434+1952 


14: 


34: 


:39 


.79 


19 


:52 


:00. 


7 


0. 


.7 


0. 


22 





.24 


0. 


25 


0. 


.24 


0. 


.6 





.9 


2 


.1 





.9 





.73 


nnyn 


+6.54c-04 


± 


6.7c-05 


2 


J1436+2321 


14: 


36: 


■AO 


.98 


23 


:21 


:03. 


2 


0. 


5 


0. 


71 





.64 


0. 


59 





56 





.4 





.8 


1 


.6 





.3 





.55 


nnyn 


+3.67c-04 


± 


3.1e-05 





J1436+6336 


14: 


36 


:45 


.80 


63 


:36 


:37. 


.8 


0. 


4 


1. 


14 


1. 


14 


1. 


.14 


1 


14 


1 


.1 





.8 





2 





.2 





.51 


nnnn 


+1.31e-04 


± 


1.2e-05 





J1437+5112 


14: 


37 


:19 


.79 


51 


:12 


:49. 


3 


0. 


.7 


0. 


.13 





.14 


0. 


15 





14 


2 


.8 


1 


.6 


1 


.7 


1 


.6 


1 


,04 


ynny 


+6.56e-04 


± 


4.6e-05 


1 


J1438+1235 


14: 


38 


:37 


.79 


12 


:35 


:34. 


2 


0. 


6 


0. 


11 





11 


0. 


10 





13 


1. 


.5 


3 


.4 


1 


.3 





.9 


1 


.27 


nynn 


+9.48e-04 


± 


1.0c+03 


3 


J1442+0625 


14: 


42 


:12 


.23 


06 


:25 


:26. 


1 


0. 


6 


0. 


10 


0. 


10 


0. 


08 





07 


3 


1 


2 


.9 


3 


.3 


4 


.4 


1 


.59 


yyyy 


+3.29e-03 


± 


5.2e-04 


2 


J1444+0257 


14: 


,14 


:31 


.76 


02 


:57 


:53. 


.4 


0. 


6 


0. 


12 





10 


0. 


11 





10 


2 


.9 


2 


.6 


2 


.9 


3 


.2 


1 


.32 


yyyy 


+2.38c-03 


± 


1.6e-04 


2 


J1444+1632 


14: 


14 


:39 


.23 


16 


:32 


:12. 


7 





.4 


0. 


10 





10 


0. 


10 


0. 


.10 


1. 


.1 


1 


.1 


2 


.1 


2 


.5 


1 


.41 


nnnn 


+4.63c-04 


± 


7.2e-05 


1 


J1446+1721 


14: 


16 


:35 


.34 


17 


:21 


:07. 


5 


0. 


3 


1. 


01 


0. 


99 





93 


0. 


81 





.6 





.2 


1 


.0 





.8 





.53 


nnnn 


+1.35c-04 


± 


2.9e-05 





J1447+1920 


14: 


:47 


:16 


.08 


19 


:20 


:46. 


9 


0. 


3 


0. 


10 


0. 


11 


0. 


11 





12 


1 


.2 


2 


.6 


3 


.1 


2 


.3 


1 


,29 


nynn 


+1.27e-03 


± 


1.5e-04 


1 


J1448+5326 


14: 


48 


:59 


.17 


53 


:26 


:09. 


2 





.7 


0. 


19 


0. 


19 


0. 


19 





19 


1 


.1 





.9 


1 


.2 





.8 





.85 


nnnn 


+6.51G-05 


± 


2.3e-05 





J1451+0106 


14: 


51 


:45 


.94 


01 


:06 


:18. 


7 


0. 


6 


0. 


08 


0. 


08 


0. 


08 


0. 


08 


2 


3 


2 


.9 


2 


.8 


1 


.4 


1 


.73 


nyyn 


+1.73c-03 


± 


1.0c+03 


3 


J1453+2648 


14: 


53 


:53 


.60 


26 


:48 


:33. 


4 


0. 


5 


0. 


35 


0. 


37 


0. 


34 


0. 


.31 





.7 


1 


.8 


1 


.8 





.8 





.63 


nyyn 


+6.61e-04 


± 


4.8e-05 


3 


J1457+0938 


14: 


:57: 


:52 


.53 


09 


:38 


:16. 


5 


1, 


1 


0. 


23 





23 


0. 


23 


0. 


23 


0. 


.4 





.8 





.8 





.3 





.75 


nnnn 


-5.92e-06 ± 1.5e-05 





J1459+2708 


14: 


59: 


:39 


.91 


27 


:08 


:18. 


8 





5 


0. 


.11 





11 


0. 


11 


0. 


.11 


1. 


.4 


2 


.6 


3 


.1 


1. 


.8 


1 


.30 


nyyn 


+9.57e-04 


± 1.3e-04 


1 


J1500+4751 


15: 


00 


:48 


.65 


47 


:51 


:15. 


5 


0. 


6 


0. 


32 


0. 


31 


0. 


33 





38 


2 


.1 


1 


.3 


1 


.4 


1 


.0 





.64 


ynyn 


+6.85c-04 


± 4.1e-05 


2 


J1502+3350 


15: 


02 


:33 


.94 


33 


:50 


:55. 


7 





4 


0. 


08 





08 


0. 


09 





10 


2. 


.1 


1 


.9 


2 


.9 


1 


.1 


1 


.55 


nnyn 


+7.01o-04 


± 1.2c-04 


1 


J1505+0326 


15: 


05 


:06 


.47 


03 


:26 


:30. 


8 


1. 


.0 


0. 


64 





81 


0. 


87 





.85 





.5 





.7 


3 


.5 





.5 





.52 


nnyn 


+1.92G-03 


± 1.6c-04 


3 


J1505+3318 


15: 


05 


:09 


.01 


33 


:18 


:32. 


5 


0. 


.4 


0. 


08 


0. 


08 


0. 


08 


0. 


.09 


1 


.8 


2 


.9 


1 


.7 


1 


.6 


1 


,61 


nnnn 


+5.11O-04 


± 7.5c-05 


1 


J1506+3730 


15: 


06 


:09 


.53 


37 


:30 


:51. 


1 


0. 


7 


0. 


76 





.70 


0. 


68 





.65 


1 


3 


2 


.9 


1 


.0 





.8 





.54 


nnnn 


+6.57c-04 


± 1.4c-04 


2 


J1514+2931 


15: 


14: 


:20 


.88 


29 


:31 


:03. 


6 


0. 


5 


0. 


12 


0. 


12 


0. 


12 





12 


1. 


.3 


1 


.7 


2 


.3 


1 


.3 


1 


.20 


nnyn 


+6.41e-04 


± 5.3e-05 


2 


J1516+1932 


15: 


16 


:56 


.79 


19 


:32 


:13. 





0. 


9 


0. 


50 





.53 


0. 


52 





49 





.9 


2 


.2 


2 


.2 


1 


.0 





.56 


nyyn 


+1.23c-03 


± 1.0e+03 


1 


J1521+1756 


15: 


21 


:17 


.58 


17 


:56 


:01. 





0. 


.8 


0. 


17 





.17 


0. 


17 





17 





.6 


1 


.1 


1 


.9 





.6 





.91 


nnyn 


+3.16c-04 


± 4.1e-05 





J1521+7858 


15: 


21 


:02 


.80 


78 


:58 


:30. 


2 





.8 


0. 


09 


0. 


09 


0. 


08 





10 


2 


.1 


2 


.0 


1 


.7 


1 


.8 


1 


.53 


nnnn 


+3.03e-04 


± 3.0e-05 






Table 4 — Continued 



Name RA(J2000) Dcc(J2000) Ha R S„(Jy) 



J1526+1632 


15: 


:26: 


41. 


90 


16: 


:32 


:45. 


9 





.8 


J1530+3758 


15 


:30: 


16. 


25 


37: 


:58: 


:31 


.1 


0, 


6 


J1534+0131 


15 


:34: 


52 


.45 


01: 


31 


01 


,2 


0. 


.8 


J1535+4836 


15 


:35: 


14. 


65 


48 


:36 


:59 


.6 





2 


J1535+6953 


15 


:35: 


19. 


22 


69: 


:53 


:18. 


3 


1 


2 


J1536+3833 


15 


:36: 


13. 


84 


38: 


33 


:28. 


,6 


0. 


5 


J1549+0237 


15 


:49: 


29. 


43 


02 


:37: 


:01 


.1 


1 


.4 


J1549+5038 


15 


:49: 


17. 


46 


50 


:38 


:05 


.7 


0. 


4 


J1550+0527 


15 


:50: 


35. 


26 


05: 


:27: 


10 


4 


1 





J1556+1230 


15 


:56: 


50. 


78 


12: 


:30 


35 


6 


1 


2 


J1559+0805 


15 


:59: 


19 


72 


08: 


:05: 


:17 


,5 


0, 


9 


J1602+3326 


16 


:02: 


07 


26 


33: 


:26 


53 





0. 


5 


J1604+1744 


16 


:04: 


26 


50 


17: 


:44 


:31 


,1 


1 





J1608+1029 


16 


:08: 


46. 


20 


10: 


:29: 


.07 


.7 


1 


.4 


J1610+7809 


16 


:10: 


50 


62 


78: 


:09 


00 


.5 


0. 


.8 


J1613+3412 


16 


:13: 


41. 


06 


34: 


:12 


:47. 


.9 





.5 


J1615+0623 


16 


:15: 


55. 


67 


06: 


:23 


19 


.0 





.8 


J1616+0459 


16 


:16: 


37. 


55 


04 


:59: 


:32 


.7 


1 


2 


J1616+4632 


16 


:16: 


03. 


76 


46: 


32 


:25. 


,2 


0. 


5 


J1617+0246 


16 


:17: 


49. 


90 


02 


:46: 


43 


.0 


1 


1 


J1619+2247 


16 


:19: 


14. 


.82 


22 


:47: 


.17 


.9 


0. 


.6 


J1623+6308 


16: 


:23: 


46. 


77 


63: 
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10.60 9.49 8.75 8.85 
0.07 0.06 0.09 0.07 
0.08 0.09 0.09 0.09 
0.10 0.10 0.10 0.10 
0.72 0.70 0.72 0.73 
0.91 0.90 0.98 0.98 
0.12 0.12 0.12 0.12 



Table 4 — Continued 



Name RA(J2000) Dcc(J2000) Ha R &,(Jy) Mod index % %err Var D(2d) tim 
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0. 
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.3 
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Name RA(J2000) Dcc(J2000) Ha R S v (3y) Mod index % %err Var D(2d) tim 
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Table 4 — Continued 



Name RA(J2000) Dcc(J2000) Ha R &,(Jy) Mod index % %err Var D(2d) tim 
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Table 5: Counts versus redshift for sources with D(2d) above the threshold 4 x 10 



Redshift 


0-0.5 


0.5-1 


1-1.5 


1.5-2.1 
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2.6-3.1 
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